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Abstract 

We introduce a new concept of quasi- Yang-Baxter algebras. The quantum quasi- 
Yang-Baxrer algebras being simple but non-trivial deformations of ordinary algebras of 
monodromy matrices realize a new type of quantum dynamical symmetries and find an 
unexpected and remarkable applications in quantum inverse scattering method (QISM). 
We show that applying to quasi- Yang-Baxter algebras the standard procedure of QISM 
one obtains new wide classes of quantum models which, being integrable (i.e. having 
enough number of commuting integrals of motion) are only quasi-exactly solvable (i.e. 
admit an algebraic Bethe ansatz solution for arbitrarily large but limited parts of the 
spectrum). These quasi-exactly solvable models naturally arise as deformations of known 
exactly solvable ones. A general theory of such deformations is proposed. The correspon- 
dence "Yangian — quasi- Yangian" and ^^XXX spin models — qnasi-XXX spin models" 
is discussed in detail. We also construct the classical conterparts of quasi- Yang-Baxter 
algebras and show that they naturally lead to new classes of classical integrable models. 
We conjecture that these models are quasi-exactly solvable in the sense of classical inverse 
scattering method, i.e. admit only partial construction of action-angle variables. 
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1 Introduction 



The quantum inverse scattering method (QISM) formulated in the works |31, 27, ^, 32, 21 
of the Leningrad group is at present time the bigest factory of quantum integrable and 
exactly solvable models most of which are of a great physical and mathematical relevance 
(see e.g. refs. [p2| , po| , |9|, 0] and references therein). The main idea of this method rests on 
the use of a special associative algebra Tr (known in the literature under names of Yang- 
Baxter algebra, Zamolodchikov algebra, Fundamental Commutation Relations, etc.)^, whose 
generators T^^, a, /3 = 1, . . . , d (forming the so-called monodromy matrix) satysfy the system 
of bilinear relations 

-Ra/375(A - fJ.)T^p{X)TsaifJ.) = T/3sifJ-)Ta-y{X)R^Spai>' - m)- (1-1) 

Tensor RafS'ysiX) (which is usually called the i?-matrix) obeys the famous Yang-Baxter equa- 
tion and completely determines the structure of algebra Tr. The role of this algebra in QISM 
can be explained as follows. By using relations (|l]^) one easily constructs the families of 
operators 

i/(A)=P[{r,,ft(A,)}] (1.2) 

commuting with each other for any values of A. Here V[...] denote some polynomials in 
generators Tq-^. (Aj) and Aj are some fixed functions of A. These operators are considered 
as integrals of motion of some quantum systems. To solve an eigenvalue problem for these 
systems one uses a special ansatz 

= Q[{Tc.Mm\o) (1-3) 

(the so-called Bethe ansatz) in which |0) is the so-called vacuum vector (playng the role 
of the lowest weight vector for the representation space of algebra Tr) and Q[...] are again 
some polynomials in generators ^'^..^^(^j) taken at certain points ^j. The coordinates of these 
points are found from a system of numerical equations known under name of Bethe ansatz 
equations. Summarizing, one can say that both the construction and solution of quantum 
problems in QISM is a purely algebraic procedure whose realization requires the knowledge 
of only two things: the defining relations of algebra Tr itself and the properties of its lowest 
weight representations. It turns out that most of the models obtained and solved in such a 
way are integrable (i.e. admit enough number of commuting integrals of motion) and exactly 
solvable (i.e. can be solved algebraically for the whole spectrum)^. This enables one to qualify 
the Tr algebra as a generator of integrable and exactly solvable quantum problems. 

In this paper we intend to introduce a new concept of "quasi- Tr algebras" . These algebras 
are distinguished by the fact that "almost" the same construction which led in the case of 
the ordinary Tr algebras to exactly solvable models, leads in the case of quasi- Tr algebras to 



Algebra Tr appeared first in Yang's papers |50[ |5l] as algebra of non-relativistic S'-matrices and in Bax- 
ter's works pi as algebra of transfer matrices. Later it reappeared in the papers by Zamolodchikov and 
Zamolodchikov [ 5^, [ssj as algebra of relativistic 5- matrices. The most regular way of introducing Tr algebras 
was formulated by Faddeev's group |, Q |l|, ||, ||, Hi and also by Drinfeld ||, |] and Jimbo [0 

who proved that Tr is a Hopf algebra. 

*In quantum mechanics the notions of integrability and exact solvability do not coincide because of the 
absence of a quantum analogue of the Liouville theorem. 
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integrable but only quasi-exactly solvable models^, i.e. models admitting an algebraic solution 
only for some limited parts of the spectrum^. The quasi- Tr algebra is a very simple but non- 
trivial deformation of the ordinary Tr algebra. The form of the commutation relations for 
Tr and quasi- Tr algebras "almost" coincide and all their properties (including construction 
of both commuting integrals of motion and their solutions) are "almost" the same. However, 
namely this "slight" difference makes the quasi- Tjj algebra a new mathematical object having 
many new exciting mathematical properties and a wide range of possible physical applications. 
The simplest example of quasi- Tr algebra has already been considered in our previous paper 



Quasi- Tr algebra is not an absolutely unexpected thing. The feeling that something like 
this may exist immediately appears if one remembers some facts from the theory of exactly 
and quasi-exactly solvable problems. Indeed, consider for example a pair of two simple 
quantum-mechanical models with hamiltonians 



H 



1=1 



dx- 



and 



^biX. 



2n 



2 2 

xf 



d 

1 - - 

2 



(1.4) 



n 



0,1,, 



(1.5) 



The first model is known to everybody. This is an ordinary d-dimensional harmonic oscillator 
which is integrable (since it admits a complete separation of variables) and exactly solvable 
(since all its eigenvalues and eigenfunctions can be constructed algebraically). The second 
model of a d-dimensional anharmonic oscillator, which appeared on the physical scene only 
recently Q| (see also [|^, ^]), is also integrable (because its variables can be separated 
in the generalized ellipsoidal coordinates p^]) but, in contrast with model (1.4), it is only 
quasi-exactly solvable (since admits only partial algebraic solution of spectral problem for any 
non-negative integer n)Q 

What can we learn from comparing these two models? First of all, we see that model 
( |1.5| ) can be considered as a deformation of the model ( |1.4D and the role of the deformation 
parameter is played by a. Second, it is clear that this deformation preserves the integrability 
property. Third, this deformation is of the splitting type, because it transforms a single 
exactly solvable model into an infinite sequence of quasi-exactly solvable models with different 



hamiltonians and different number of exact solutions [46|. 



It turns out that this situation is quite general. In fact, it can be shown that any inte- 
grable and exactly solvable model satisfying some simple and rather general conditions can 
be constructively deformed into an infinite sequence of integrable and quasi-exactly solvable 
models 14^. In particular, this is true for models associated with various Tr algebras and 



^The notion of "quasi-exact solvability" was introduced in paper pq ]. 

^Numerous examples of quasi-exactly solvable models and various methods for their construction can be 
found in several review articles |6|, |38| ^] and in the book . 

^Indeed, it can be explicitly shown that, for any given n, the number of algebraically constructable solutions 
of model O) is only (n + d)!/(n!d!). 
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obtainable by means of QISMg. But if so, may be this deformation can be understood on a 
purely algebraic level? May be, there exists some non-trivial deformation of the underlying 
Tr algebra, which could be called "quasi-T^ algebra", and which would naturally lead to 
quasi-exactly solvable versions of integrable models? 

The aim of this paper is to give a positive answer to this question. To make the text more 
readable, we devote the main attention to the simplest rational T^-algebra corresponding 
to a monodromy matrix of the size 2x2. The rationality means that both generators 
Tafii^), a, (3 = 1,2 and tensor Rap^ysW, a, /?, 7, (5 = 1, 2 are rational functions of A. Such an 
algebra is called the Yangian 3^[s/(2)] and is characterized by the following choice of non-zero 
elements of tensor Raf3ysW 

-Rllll(A) = i?2222(A) = 1, 

i?i2i2(A) = i?2i2i(A) = A/(A + r?), 

i2i22i(A) = i?2ii2(A) = r?/(A + r/), (1.6) 

where is a parameter. For 77 7^ the Yangian 3^[s/(2)] is a non-commutative and non- 
cocommutative Hopf algebra H], ^ . In the limit r/ — > it degenerates into a special (infinite- 



dimensional) Lie - Gaudin algebra |13] which we denote by ^[s/(2)]. The completely inte- 
grable and exactly solvable models associated with algebras ^[sZ(2)] and 3^[s/(2)] are respec- 
tively known under names of Gaudin and XXX spin models. 

The paper is organized as follows. After short exposition of properties of ordinary algebras 
^[sZ(2)] and 3^[s/(2)] (sections ^ and ^) we explain the reader what we mean under quasi- 
t/[sZ(2)] and quasi-3^[s/(2)] algebras and present their defining relations (sections |^ and P) . We 
also give a proof of the existence of these algebras by contructing their explicit realizations 
(sections |^ and 0). In particular, we demonstrate that they actually can be considered as 
deformations of algebras ^[sZ(2)] and 3^[s/(2)] and show that these deformations are of the 
splitting type, i.e. the generators of "quasi" algebras contain an additional subscript n. The 
non-triviality lies however in the way in which this subscript appears in the defining relations 



(1.1) which, as for the rest, are the same as in the undeformed case. We show also that 
the proposed deformation preserves the integrability property (sections ^ and ^). Following 
general prescriptions of QISM, we construct the infinite series of the associated integrable and 
quasi-exactly solvable models (sections and ^jl^) and present their partial solutions in the 
framework of the algebraic Bethe ansatz (sections 4 and 10). We call the models obtained in 
such a way the quasi-Gaudin and quasi-XXX models, respectively. 



In sections |13| - |16| we discuss the obtained results from the point of view of the existing 
methods in the theory of quasi-exactly solvable systems. We consider three general methods 
known under names of "partial algebraization method" (section |l^, "inverse method of 
separation of variables" (section |l^) and "projection method" (section |l^) , and show that, in 
contrast with the first and second methods, which are either not applicable or only partially 
applicable to the study of quasi- Tr algebras, the third one enables one to reproduce correctly 
practically all features of the proposed formalizm. However, the mathematical technique used 
in this method is so cumbersome and model-dependent that it hardly can be practically used 
beyond the simplest cases of algebras t/[s/(2)] and 3^[s/(2)] and their "quasi" analogues. 

*The case of the simplest models of such a sort (the so-called Gaudin models |l^) was considered in refs. 
[ ^sf |46| . For more general discussion of this subject see sections |l^ and |l6| of present paper. 
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Fortunately, the "quasi-QISM" formalism which we propose in this paper is free of all 
these deficiences. Indeed, in sections |l^ - ^ we show that this formalism is quite general and 
applicable not only to Q[sl{2)] and 3^[s/(2)] algebras, but also to all Tr algebras, irrespective 
of their dimension and a concrete form of tensor Ra/3^sW- For this we will use a model 
independent formalizm based on the notion of Z'^-graded unital associative algebras (which 
in this paper we call simply A-algebras) and construct quasi-deformations immediately for the 
latter. The obtained results are automatically applicable to any Tr algebras since they can 
be considered as particular cases of ^-algebras. This gives us the hope that using this general 
deformation technique we will find a lot of new quasi-exactly solvable models interesting both 
from physical and mathematical points of view. 

One of questions discussed in this paper (see sections ^ and is related to a special 
limit of quasi- Tfj algebras when the so-called "spectral parameter" A characterizing their 
generators tends to infinity. We know that in the case of ordinary (classical or quantum) Tr 
algebras such a limit leads respectively to ordinary finite-dimensional Lie algebras or their 
(7-deformed (Licq) versions. It turns out that performing the same limiting procedure with 
quasi- Tfj algebras we also obtain some new interesting mathematical structures which can be 
called quasi-Lie and quasi-Lieg algebras. 

In section ^ we introduce the concept of classical quasi- Tr algebras. These algebras 
naturally arise as the result of dequantization of quantum quasi- Tr algebras discussed in 
preceding sections. Their generators are realized as functions on the phase space and obey 
the special classical commutation relations. The latter are expressed through a special non- 
standard (quasi-Poisson) bracket which can be considered as a deformation of the ordinary 
Poisson bracket. We show that the the role of classical quasi- Tr algebras in the classical 
inverse scattering method (CISM) is the same as that of quantum quasi-T^ algebras in QISM. 
The classical models associated with classical quasi-T^ algebras are integrable but it seems 
that they are only partially solvable in the framework of CISM. 



In the concluding section 23 we give a list of problems which naturally arise from the 



results of the present paper and, in our opinion, deserve a very careful study. 



2 Gaudin algebra and Gaudin model 



In this section we remind the reader some basic facts concerning Gaudin algebras, their 
representations and properties of the associated Gaudin models. More detailed exposition of 
this subject can be found in refs. ||, ||, ||, |6|, 

As we mentioned in Introduction, the Gaudin algebra C/[s/(2)] can be obtained from the 
Yangian y[sl{2)] (which is characterized by i?-matrix of the form (|1.6D and commutation 
relations (0) ) in the limit when ry ^ 0. For this the generators TapiX), a,/3 = 1,2 should 
be taken in the form 



rn(A) = l + r?50(A) + O(r?2), 

T22{X) = l-r?50(A) + O(r?2), 

ri2(A) = -rjS-{X) + 0{v'), 

rn(A) = +^5+(A) + 0(r,2), 



(2.1) 



6 



where S^{X), S~{X) and ^^(A) are just the generators of ^[s/(2)]. These generators, which are 
respectively called neutral, lowering and raising generators, obey the following commutation 
relations 

5°(A)5°(^) -5°(/i)5°(A) = 0, 
5+(A)5+(;u) -5+(/x)5+(A) = 0, 
S-{X)S-{f,)-S-ifi)S-{X) = 0, 

S^{X)S+{fi)-S+{fi)S\X) = -{X- n)-^{S+{X) - S+{n)}, 
S'{X)S-if,)-S-{f,)S\X) = +{X-f,)-'{S-{X)-S-if,)}, 
5-(A)5+(/.)-5+(/x)5-(A) = -2(A-/.)-H5°(A)-50(/.)}. (2.2) 

From ( ^.21) it immediately follows that operators 

H{X) = S\X)S\X) - ^5-(A)5+(A) - ^5+(A)5-(A) (2.3) 
form a commutative family, 

[HiX),Hifi)]=0, (2.4) 

and thus can be interpreted as integrals of motion of a certain quantum completely integrable 
model^. The latter is known under name of Gaudin model. 

The role of the carrier space in which the operators H{X) actP| is played by the represen- 
tation space of Gaudin algebra. In order to construct it one needs to fix the lowest weight 
vector |0) and the lowest weight function F(X) obeying the relations 

50(A)|0)=F(A)|0), 5-(A)|0)=0. (2.5) 

After this, we can define the representation space as a linear hull of vectors 

\^l,...,Cn) = S+{Cn)...S+iCl)\0), (2.6) 

with arbitrary n and ^i, . . . , We shall denote this space by Wp(^xy 
The "Schrodinger equation" for the Gaudin model reads now 

H{X)(t> = E{X)4>, ^eWp^xy (2.7) 

The beauty of this equation lies in the fact that all its solutions can be obtained algebraically 
within the so-called Bethe ansatz 

(/> = 5+(go5+ (e„„i) . . . 5+(e2)5+(ei)io), (2.8) 

^Using concrete realizations of generators of algebra G[sl{2)], one can easily check that the number of 
independent operators belonging to the class H{X) is sufficiently large for claiming that this model is completely 
integrable. 

^''We try to avoid the use of the term "Hilbert space" because we do not intend to discuss in this paper the 
hermiticity properties of operators H{X) and the normalizability of their eigenstates. 
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in which n is an arbitrarily fixed non- negative integer and ■ ■ ■ ,Cn are some still unknown 
numbers. To demonstrate this, one needs to act on Bethe vector by operator H{X) and, using 
commutation relations 

H{\)S+{Ci) = S+{Ci)H{X) + 2 ^^fe)'g°(^)-f^(^)^°fe) (2.9) 

which elementary follow from ( |2.2| ), transfer H{X) to the right. The neutral operators S^{X) 
and S'^{S,i) appearing in the right hand side of ( |2.9D also should be transfered to the right. 
For this the fourth relation in the list ( |2.2[ ) can be used. Finally, we obtain a number of terms 
with operators H{X), S^{X) and S^{(,i), i = 1, . . . ,n acting immediately on the lowest weight 
vector |0). After this one should get rid of these operators by using formulas ( |2.5[ ) together 
with the relation 

H{X)\0) = {FHX)+F'{Xm (2.10) 

which elementary follows from (|2.5| ) and (|2.3DP^ . The result will consist of two parts. The first 
part will be proportional to the Bethe vector while the second part will contain n terms 
not having the initial Bethe form (these are the so-called "unwanted terms" ) . The condition 
of cancellation of these terms imposes some special conditions on numbers ^j, i = 1, . . . ,n 
which up to now were considered as free parameters. An explicit form of these conditions, 
which are commonly known as Bethe ansatz equations, reads 

n ^ 

E F3^ + ^te) = 0' i = l,...,n. (2.11) 

Putting all these facts together, one can easily obtain the final solution of the Gaudin spectral 
problem: 

E{X)=F\X)+F'{X) + 2J2^^^^T^P^- (2.12) 

i=i « 

It is obvious that if F{X) is a rational function, then the number of solutions of Bethe 
ansatz equations ( 2. 11]) is finite for any n. But since n is an arbitrary non-negative integer. 



the complete number of solutions of problem (|2.7| ) is infinite. It can be shown that these 
solutions expire all possible solutions of problem ( p. 7] ) and therefore it can be qualified as 
exactly solvable. 



3 Quasi- Gaudin algebra and quasi- Gaudin models 



Let us now try to explain the reader what we mean under words "quasi-Gaudin algebra" . Note 
that this will be namely an attempt of an explanation and not a list of rigorous definitions 
and theorems. A strict mathematical interpretation of formulas of this section will be given 



later in sections 17 



^^In order to derive ( |2.1C| ) one should first take the last relation in the list (2.2) in the limit when ^ ^ A, 
and using it rewrite the operator (Ei) in the form H{X) = [S°{\)f + [5° (AY]' - S'+(A)5"(A). After this one 



should act by this operator on the vacuum vector and then use formulas (2.5) 
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Let S^{X), (A) and S^{X) denote the operators parametrized by an integer number n 
and obeying the relations 

50(A)50(^)-50(^)50(A) = 0, 

5',7-i(A)S'~(/i) - S'^_i(m)5'^(A) = 0, 

50_,,(A)5+(m)-S+(/.)50(A) = -(A-^)-HS+(A)-S+(^)}, 

50_i(A)5-(m)-S-(^)50(A) = +(A-^)-HS-(A)-S-(^)}, 

S-+i(A)S+(^)-5+_i(m)S-(A) = -2(A-^)-iK(A)-S0(/.)}. (3.1) 



We consider (3.1) as the definig relations for a certain infinite-dimensional algebra which is 
very similar to the Gaudin algebra but, in contrast with the latter, it is not a Lie algebra. 
We shall call it "quasi-^[s/(2)] algebra" and the expressions staying in the left hand sides 
of relations ( |3.1| ) we refer to as "quasi-commutators" . A most unusual feature of algebra 
Q[sl{2)] is that it is partially free i.e. does not contain quasi-commutation relations between 
generators whose indices differ more than by one. One can also say that this alhebra is local 
in the space of a discrete parameter n. Note that the form of quasi-commutation relations is 
invariant under special "quasi-similarity" transformations 

5°(A) ^ UnSl{X)U-\ 

St{X) ^ Un+lStAmn\ (3.2) 

in which [/„ are arbitrary invertible operators. 

One of the remarkable properties of quasi-^[sZ(2)] algebra is that its generators (exactly 
as in the ^[s/(2)] algebraic case) can be used for building hamiltonians of some completely 
integrable quantum systems. Indeed, using quasi-commutation relations (^]^) it is not difficult 
to prove that the operators 

H^{\) = 50(A)50(A) - \s-^^{\)St{\) - \st^^{\)S-{\) (3.3) 
form commutative families for any n: 

[Hn{X),Hn{tl)]=0. (3.4) 

At the same time, the operators Hn{X) and HmifJ-) with different n and m do not generally 
commute with each other because of the absence of appropriate (non-local) relations in algebra 
Q[sl{2)]. This means that formula (|3.3| ) describes in fact an infinite set of different integrable 
model!:^. This is, may be, the most important difference with the Gaudin case where an 
analogous construction leads to a single completely integrable model — the Gaudin model. 

Exactly as in the case of Gaudin algebra, the final construction of models with hamil- 
tonians (|3.3|) needs the specification of a carrier space in which these hamiltonians act. It 



^^Later, in section after constructing a concrete realization of algebra C7[s^(2)] we will be able to demon- 
strate that operators H„{X) and Hm{fJ,) actually do not commute with each other. At the same time, the 
number of independent operators contained in the family i/,i(A) with a given n is sufficiently large to claim 
that the corresponding models are completely integrable. 
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would be natural to identify it with the "representation space" of our algebra]^. The latter 
can be constructed almost in the same way as in the Gaudin case. There is however an 
essential difference which lies in the fact that the quasi-analogoues of relations ( |2.5D cannot 
be written simultaneously for all "neutral" and "lowering" operators 5*^ (A) and S~{X), i.e. 
for any n. This simply would lead to an overdetermined system of equations for the vacuum 
vector |0). Roughly speaking, to have the same number of equations for |0) as in (|2.5D , one 
should restrict ourselves to one arbitrarily fixed value of n. Taking for example n = 0, we 
obtain 

50(A)|0)=F(A)|0), 5-(A)|0)=0, (3.5) 

where -F(A) is, as before, a certain function playing the role of the lowest weight. The 
imposibility of writing the relations like (^) for any n is another unusual property of quasi- 
^[sZ(2)] algebra, and it, strictly speaking, should be postulated from very beginnin^^. The 
presence of such a postulate will play a very important role in our further considerations. 

After fixing the "lowest weight vector" |0) together with the corresponding "lowest weight 



function" F{X), we can define the "representation space" by analogy with (2.6) i.e. as a linear 
hull of vectors obtained after action of operators S^{^i) on |0). It is however intuitivelly 
clear that there is a big freedom in choosing such vectors because the indices of raising 
operators may take infinitely many different values. In order to restrict this freedom (which, 
unfortunately, we cannot control at present time), we need a "guiding principle", which, in 
particular, could be based on a selectional use of some symmetries. Looking for example 
at the vectors (|2.6| ) forming the representation space of ordinary Gaudin algebra, we see 
that, due to the commutativity of raising generators 5'''(^i), . . . , 5''^(^„), they are symmetric 
under all permutations of numbers ^i, . . . It is natural to demand the same symmetry 
in the case of quasi-Gaudin algebras. Now, however, the raising operators do not commute 
anymore, they only "quasi-commute" , and therefore the only way of having such a symmetry 
is to choose the corresponding vectors in the form 

16, • • ■,^m)k = S:^^i.{^rn)S:^_i_^_k{^m^l) • • • S'^+fe (?2)5'jt^A: ) l^) (^-S) 

with arbitrary k £ Z , m £ N and 6, . . . , As before, we denote the "representation space" 
of algebra ( |3.1D defined in such a way by Wp(^xy 



Now we can complete the construction of integrable models (3^) associated with algebra 



(|3.1|). We postulate that "Schrodinger equations" for these models read 

i7„(A)0 = S(A)(/>, (/>gVFf(A), n = 0,1,2,.... (3.7) 

Hereafter we shall call such models the "quasi-Gaudin models" . 

^''Generally speaking, the object which we intend to call the "representation space" is not a true repre- 
sentation space for quasi-C/[sZ(2)] algebra since, even after choosing an appropriate basis in it, we cannot 
represent any element of this algebra in a matrix form. This is possible only for some special realizations of 
quasi-C/[sZ(2)], one of which will be considered in section ^. 

^''Later, in section ^ we demonstrate that the lowest weight vector |0) is actually not annihilated by the 
"lowering operators" 5,7 (A) if n 7^ 0. Here, however, we suggest the reader to consider this fact as a postulate. 
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4 Bethe ansatz for quasi-Gaudin models 



It turns out that, in full analogy with the Gaudin case, the solutions of spectral equations (3.7) 



can be found by means of the algebraic Bethe ansatz. But before constructing the appropriate 
Bethe vectors one should discuss some auxiliary problems concerning the properties of quasi- 
Gaudin algebra. 

It is known that in the Gaudin case the lowest weight vector is always a solution of the 
Gaudin spectral problem (see for example formula ( ^.lOD ). Is this true in the case of quasi- 
Gaudin models? In order to answer this question, one should simply act by the operator 
Hn{X) on the vector |0) and look at the result. For this it is convenient to rewrite the 
operator Hn{X) in a little bit different form. Using the last formula in the list ( fj.l[ ) and 
taking in it the limit /u — > A we obtain 

5-+i(A)5+(A) = S+_,{X)S-{X) - 2|^50(A). (4.1) 



Substituting (4.1) into expression (|3.3|) for Hn{X), we get the needed expression for Hn{X) 



Hn{X) = S0(A)50(A) + ^S^^iX) - S+_,{X)S~{X). (4.2) 

Remembering now that for n ^ the vacuum vector |0) is not in general annihilated by 
the lowering operators ^"(A) and is not an eigenvector of neutral operators ^^(A), we can 
conclude that it may be an eigenvector of (|4.2| ) only if n = 0. Using formulas (|3.5D , we obtain 
for this case the relation 

Ho{X)\0) = {f\X) + F'{X))\0) (4.3) 
which is similar to the Gaudin relation ( 2.1C| ). 



Let us now return to the main problem of solving equations ( |3.7D by means of Bethe 
ansatz. It is natural to take the Bethe vector in the form 

^ = s++,(e™)5+_i+fe(e,n-i) . . . ^2^(6)5,^1(6)10) (4.4) 

with some k Z and m N, and, using the relations 

U f\\C+ (t \ C+ (C \U (\\ I n,^n-\^in)S%^-\{X^ — S^-\{X)S^_^{(,r^ 

Hn{X)S^_^{^n) = 5^_i(C„)F„_.i(A) + 2 — , (4.5) 

which follow from the basic relations (p.l|), try to transfer the operator Hn{X) to the right 
(exactly as in the ordinary Gaudin case). From formula (4.5) it is seen that, in order to start 
performing this procedure, it is necessary to take 

m + k = n — 1. (4.6) 

Otherwise we would be unable to transfer the operator Hn{X) to the right because of the 
absence of the corresponding commutation relations in ( |4.5| ). From formula ( |4.5| ) it is readily 
seen that permutation of the operator Hn{X) with the first raising generator S^_i{£^n) of 
Bethe vector replaces Hn{X) by Hn-i{X). After this, Hn-i{X) will appear in front of the next 
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raising generator S^_2{Cn-i), and we can again permute these operators by means of the same 



formula (4.5) but with n replaced by n— 1. Obviously, this procedure can be continued further. 



The same relates to the neutral generators ^^.^(A) and S^_i{^n) appearing in the right hand 



side of (4^) and thus acting immediately on S^_2{^n-i)- The fourth relation in the list (3.1) 
enables us to perform the permutation of these operators, after which the neutral generators 
appear in front of the generator S^_^{^n-2) and will have the index n — 2. This enables one 
to use again the forth formula in (|3.lD , etc. Summarising, one can say that any permutation 
of H- and S'^'-operators with raising generators forming the Bethe vector decreases their 
indices by one. Finally, when all the transferences will be completed, we obtain a number 
of terms with operators Hn-mW, S^_^{X) and S^_^{^i), i = I, . . . ,n acting immediately 
on the lowest weight vector |0). The standard presciptions to Bethe ansatz technique (see 
e.g. section imply that all these operators should be absorbed by this vector. But such 
an absorbtion is possible only if the lowest weight vector is an eigenvector of these operators, 
which, according to the previous results, is possible only if 

n - m = 0. (4.7) 

Comparing formulas ( [4.6D and ( ^.7] ) we can conclude that the only case when the ansatz ( [4.4] ) 



for equation (|3.7| ) may lead to some algebraic solutions corresponds to the choice 

k = —1, m = n. (4-8) 
After using the restrictions (^]^) the ansatz (|4.4| ) takes the following final form 

<A = 5+_i(^„)5+_2(en-i) . . . St{i2)SUiim- (4.9) 

Now it remains only to check that it actually contains solutions of the quasi-Gaudin spectral 
problem (|3.7| ). This can be demonstrated exactly in the same way as in the ordinary Gaudin 
case. Getting rid of the operators Hq{X), Sq{X) and 5o(^j), i = 1, . . . ,n standing in front of 
the vacuum vector by means of formulas (|3.5| ) and ( [4.3]) , we obtain two sorts of terms. One 
of these terms will be proportional to the Bethe vector (j), while the other ones will represent 
the good old "unwanted terms" . It is not difficult to show that the condition of cancellation 
of the latter reads 

n 

E F3^ + ^te)=0' i = l,...,n, (4.10) 
and the final expression for the eigenvalues E{X) is given by 

E{X) = F\X) + F'(A) + 2 V (4.11) 



This gives us the desired algebraic solution of problem (|3.7| ). From the above consideration 
it is clearly seen that models with hamiltonians Hn{X) are typical quasi-exactly solvable 
ones. Indeed, the linear hull of all admissible Bethe vectors (|4.9| ) forms only a certain small 
subspace in the space of vectors of the general form ([4.4|). The number of explicit solutions 
of problems (^3) is finite for any given n if F{X) is a rational function. At the same time 
the dimension of the space Wp(^x^ is generally infinite-dimensional. It is interesting that we 
obtained exactly the same set of solutions as in the case of ordinary Gaudin model. But now 
these solutions are distributed between different models. 
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5 The existence of quasi-Gaudin algebras 



Everything what we said in sections ^ and Q can be checked by means of direct calculations. 
The only thing which may seem dubious is a very existence of an algebra with such strange 
properties. In order to make sure that it does actually exist, one should construct at least one 
of its explicit realizations. One of possible realization can be constructed immediately from 
the generators of Gaudin algebra. Let us take the generators S^{X), S~{X), 5"+ (A) introduced 
in section |2| and assume that they act in the representation space characterized by the lowest 
weight function F^{\). Let us also introduce a special limiting operator 

S° = lim A5°(A). (5.1) 

A^oo 

Then the desired realization reads 

F" - 5° + n 
A — a ' 
FO - S° + n + 6 



5- (A) 


= s-{x) + 




= 50(A) + 


S-{X) 


= S~{\) + 



A — a 
FO - 5° + n + 26 



(5.2) 



A — a 

where a and b are arbitrary complex parameters and = limA^oo AF'^(A). It can be easily 



checked that operators (5^) actually satisfy the quasi-commutation relations (3.1) 



The representation space of this algebra coincides by construction with the representation 
space of the Gaudin algebra. The lowest weight vector |0) is the same as in section ^ which 
gives us the possibility of checking formulas ( ^.5]) by using their Gaudin analogues (|2.5D . 
Taking into account that 

5°|0)=F°|0) (5.3) 

and using formulas ( |2.5| ), we obtain: 

5°(A)|0) = (fO(A) + ^) |0), S~iXm = |0). (5.4) 

We see that the lowest weight vector is actually not annihilated by "lowering" operators if 
n 7^ 0. Taking n = in ( p. 41 ) we obtain the expression for the lowest weight of quasi-Gaudin 
algebra 

F(A)=fO(A) + -^. (5.5) 
A — a 



Now it is absolutely clear that algebra ( |3.1D is a continuous deformation of the Gaudin algebra. 
The role of the deformation parameter is played by a. If a — > oo, then the n-dependence 
of the operators ( |5.2| ) disappears and they transform into the ordinary generators of Gaudin 
algebra. Respectively, the quasi-commutators in ( |3.lD become the ordinary ones. As to the 
formulas ( |3.5D defining the representations of algebra ( |3.1| ), they, in the limit a — > oo, also 
transform into relations (|2.5[) for the Gaudin algebra. 
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6 Quasi-exact solvability of quasi- Gaudin models 

Substituting ( ^.2D into ( |3.3D , we obtain the following explicit form of quasi-Gaudin hamilto- 
nians: 

Hn{X) = S\\)S\\)-\s-{\)S+{\)-\s+{\)S~{\) + 

25°(A)(n + fe + F° - 5°) - S-{\){n + 2b + F^- 5°) - S+{X){n + F° - g°) 

A — a 



Now we can check in an independent way that the operators ( |6.1[ ) actually describe quasi- 
exactly solvable models. Denote by a linear hull of all vectors S~^{S,k) ■ ■ ■ 'S''''(^i)|0) with 
arbitrary ■ ■ ■ ,Ck aiid k < n. It is known that dim < oo for any n if F^{X) is a rational 
function. From the obvious relations S^{X)^n C <I>„±i, ^"(A)^^ C ^n, {n + - S^)^n C 
<I>„_i and {n-\-F^ — S^)^m C for m 7^ n it immediately follows that each of the operators 
Hn{\) admits only one invariant subspace, This subspace is finite-dimensional and 



therefore the models (|6.l| ) are quasi-exactly solvable in the standard sense of this word (for 



the "standard" quasi-exact solvability see e.g. refs. |2^, p6[). 

7 Quasi-s/(2) algebra 

It is known that the sl{2) algebra can be considered as a limiting case of the Gaudin ^[s/(2)] 
algebra ( |2.1| ). Its generators , and are given by 

5° = lim AS°(A), = lim \S-{\), S+ = lim A5+(A), (7.1) 

A— >cxD X^oo A— >oo 

and satisfy the usual commutation relations 

s^s--s-s^ = -S-, 

S-S+-S^S- = 25°. (7.2) 

It is natural to ask what kind of algebra will appear if we consider an analogous limit of 
quasi-t/[s/(2)] algebra? By analogy with the Gaudin case, we can define the generators of 
this algebra by the formulas 

5° = hm A5°(A), S- = lim XS~{X), S+ = lim A5+(A). (7.3) 

A— >oo A-^00 A^oo 



Multiplying the quasi-commutation relations (3J.) by A^ and tending both A and fi to infinity 



we easily derive the relations between these generators, which read 



CU c— 0— cU c— 

'-'n— I'-'n *-'n '-'n '-'n ; 

qO c+ 0+ qO I Q+ 

^n+l^n ~ ^n-l^n = 25°. (7.4) 
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We can consider (|7.4D as the definig relations of a certain modification of the sl{2) algebra 
which can be called "quasi-s^(2) algebra". Despite the fact that it is not a Lie algebra, it has 
many properties similar to those of the ordinary sl{2). For example, it has a quasi-analogue 
of the Casimir operator, 

which quasi-commutes with all generators: 

CnS^_l = S^_iCn-l, CnS^^i = 5'„_|_]^Cn+l, CnS^ = S^Cn- (7-6) 

The "representations" of quasi-s/(2) algebra can be constructed in the same way as in the 
quasi-Gaudin case. Defining the "lowest weight vector" |0) by 

5o»=F|0), S,\0) = 0, (7.7) 

where is a certain arbitrarily fixed number, we can define the "representation space" Wp 
as a linear hull of vectors 

\in)k = 5'^+fc5'^_i+fc • • • S2+kSt+k\^), (7.8) 

with arbitrary k Z and m ^ N . If we are interested in the spectrum of quasi-Casimir 
operator, then it is easily seen that, in contrast with the standard sl{2) case, it is not infinitely- 
degenerate and contains only one exactly constructable eigenvector. For example, the only 
possibility for vector ( [7. 81) to be an eigenvector of (|7.5| ) is realized when k = —1 and m = n. 
In this case, the corresponding eigenvalue is equal to F{F — 1). So, one can say that the 
quasi-Casimir operators represent the simplest (and, in some sense, trivial) quasi-exactly 
solvable models. 

It is remarkable that the generators of quasi-sZ(2) algebra can be realized in terms of 
the generators of ordinary sl{2) algebra. The corresponding formulas can be obtained after 
substituting formulas (|5.2|) into (|7.3|) and read 



= 5^ + n + F°-5°, 5° = n + 6 + F°, 5+ = 5+ + n + 26 + F° - 5°. (7.9) 
This realization corresponds to the choice F = -\-h \n ([7.7]). 

8 The Yangian and the XXX model 

Let us give some brief exposition of facts concerning the Yangian 3^[s/(2)], its representations 
and properties of the so-called XXX models associated with it. A more detailed exposition 
of this subject can be found in refs. 2C]. The commutation relations for generators 



2o;/3(A), a, /3 = 1,2 can be extracted from formulas ( |1.1[ ) and ( |1.6D . Introducing the commonly 
used notations 

rn(A) = A(A), ri2(A) = i?(A), r2i(A)=C(A), r22(A)=D(A), (8.1) 
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we can rewrite these relations in more explicit form 



(A 


- f,){A{X)B{i,) 


-i?(/x)A(A)} 


= v{A{fJ.)B{X) 


- A{X)Bifi)} 


(A 


- f,){A{X)C{^l) 


-C(/.)A(A)} 


= v{A{X)C{f,) 


- Mf^)CiX)} 


(A- 


- fi){D{X)B{i^) 


-Bifi)DiX)} 


= v{D{X)B{fi) 


-Di^i)BiX)} 


(A- 


- ^l){D{x)ci^,) 


-Cifi)D{X)} 


= v{DifJ.)C{X) 


- DiX)Cif,)} 


(A- 


- l,){A{X)Dit,) 


-Dif,)AiX)} 


= v{C{fJ.)B{X) 


-CiX)B{f,)} 


(A- 


- ^l){B{x)Cif,) 


-C{f,)DiX)} 


= v{Difi)AiX) 


-DiX)Aif,)}. 



Using these relations one can easily prove that the operators 

H{X) = A{X) + D{X) 

commute with each other, 

[HiX),H{f,)]=0, 



(8.3) 



(8.4) 



for any A and /x. This enables one to consider them as integrals of motion of a completely 
integrable quantum system which is known under name of XXX mode|^. The hamiltonians 
H{X) of these models act in the representation space of algebra 3^[s/(2)] which is defined as 
a linear hull of vectors 

\^i,...,Cn) = C{U---C{Cim (8.5) 

with arbitrary n and ^i, . . . , Here |0) is the lowest weight vector, satisfying the conditions 

i?(A)|0)=0, .4(A)|0) =a(A)|0), I)(A)|0) = d(A)|0), (8.6) 

and functions a(A) and b(X) play the role of lowest weight functions determining the form of 
the represent ationf^ Denoting the representation space defined by formulas ( |8.5| ) and (|8.6| ) 
by W^a(A),d{A)5 we can write down the spectral problem for ^^(A) 

H{X)<l> = E{X)(t>, 0e W,(A),rf(A). (8.7) 

It turns out that this problem is exactly solvable within the Bethe ansatz 

</. = C7(C„)...C(6)lO). (8.8) 

This can be demonstrated exactly as in the Gaudin case. We simply should act on ( ^.8]) by 
operator H{X) and transfer the operators ^(A) and D{X) (from which it consists) to the right 
by using bilinear relations 



AiX)C{^i 
D{X)CiC,) 



A-^i 



^ c(A)A(e.) + '-^^^^ 



A-e. 



C{X)D{i.,) + 



A-e* 



CmA{X), 
C{ii)D{X), 



^^More precisely — the XXX inhomogeneous magnetic chain. 

^^Some formulas of this section differ a little bit from those used in the standerd literature since, instead 
of the standardly used highest weight representation we are considering here the lowest weight one. But the 
results are, of course, the same. 
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which triviahy follow from (8^). As soon as operators A{X), D{X) (as well as operators 
j4(^j), Z)(^j), i = l,...,n arising after these transferences) appear in front of the lowest 
weight vector |0), they should be absorbed by it. The Bethe ansatz equations guaranteeing 
the absence of the unwanted terms have the form 

and the final expression for the spectrum of the XXX model reads 

m = a{X) n + d{X) n (8.11) 

i=l ^» i=l ^ 

If, for example, a(A) and d{X) are rational functions, then the number of solutions of Bethe 
ansatz equations ( ^.10[ ) is finite for any n. But since n is an arbitrary non-negative integer, 
the complete number of solutions of problem ( |8.7|) i s infinite. It can be shown that these 
solutions expire all possible solutions of problem ( |8.7D and therefore it is exactly solvable. 

9 Quasi- Yangian and quasi-XXX models 

Let us now construct the quasi- Yangian. Its generators depend additionally on an integer n 
and can be denoted by An{X), i?„(A), C„(A) and Dn{X). The quasi-commutation relations 
for these generators read 

(A - f,){An-liX)Bnifi) - Bnifi)AniX)} = (^)B„(A) - (A)5„ (;u) } 

(A - fi){An+l{X)Cn{fl) - Cn{^i)An{X)} = r/{^„+i (A)C„(^) - An+l{fl)Cn{X)} 

(A - fi){Dn-i{X)Bn{fi) - BMDniX)} = r?{D„„i (A)S„ (^) - Z)„_i(^)5„(A)} 

(A - fl){Dn+l{X)Cn{fl) - Cn{^l)Dn{X)} = (//)C„(A) - Dn+i{X)CM} 

(A - fi){An{X)Dn{fi) - Dn{^l)An{X)} = r?{C„_i (//)S„ (A) - C„_i(A)S„(/i)} 

(A - fi){Bn+l{X)Cn{fl) - Cn-l{fl)Bn{X)} = v{Dn{f^)AnW " Dn{X)AM}. (9.1) 

It is not difficult to check that these relations are invariant under "quasi-similarity" transfor- 
mations 

^„(A) ^ UnAn{X)U-\ 
Bn{X) ^ Un-lBn{X)U-\ 
CnW Un+lCn{X)Un^ , 

Dn{X) ^ UnDn{X)U~\ (9.2) 

in which C/„ are arbitrary invertible operators. 

Using ( p. ID one can easily prove that the operators 

Hn{X) = An{X) + Dn{X) (9.3) 

commute with each other, 

[Hn{X),HM]=0, (9.4) 
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for any A and /j, and any given n. However, the -ff-operators with different indices n do 
not generally commute with each other. This leads us to a conclusion that we deal with 
an infinite series of different completely integrable quantum systems which we shall call the 
quasi-XXX-models. 

In order to formulate a spectral problem for hamiltonians Hn{\) one should construct a 
carrier space in which they act. This space, as usually, can be identified with a "representation 
space" of quasi-3^[sZ(2)] algebra. For constructing the latter we can use the same reasonings 
as in section ^ and define it as a linear hull of vectors 

I'^l, • • • ■,im)k = Crn+k{im)CTn+k-l{im-l) ■ ■ ■ C'fc+2 (C2)C'fc+l ('^l) |0) (9.5) 

with arbitrary m ^ N , k €z Z and Cii • • • > Cm- It is not difficult to see that, due to commutation 
relations ( |9.lD , the vectors (|9.5| ) are symmetric with respect to all permutations of numbers 
■^1) • • • ! Cm- Here, as before, |0) is the lowest weight vector, satisfying the restricted conditions 

BoiXm=0, Ao(A)|0) =a(A)|0), Z)o(A)|0) = d(A)|0), (9.6) 

and functions a(A) and d{X) play the role of the lowest weight functions determining the form 
of the "representation". We denote such a "representation space" by Wa(x),d{x)- In the next 
section we demonstrate that the lowest weight vector is actually not annihilated by other 
lowering operators Bn{X) with n ^ 0. Moreover, we show that it is even not an eigenvector 
of operators Dn{X) if n 7^ 0. In conclusion we present the final form of the spectral problem 
for operators Hn{X), which reads 

Hn{X)<p = E{X)cl>, <p G W^,(A),rf(A) (9.7) 

and whose solutions will be discussed in the next section. 



10 Bethe ansatz for quasi-XXX model 

The results of section |^ enable one to show that spectral problems (|9.7| ) can be algebraically 
solved within the Bethe ansatz 

'P — CjTt+fc ('^m)Cm+fc— 1 ('^m— 1 )---Ck+2ia)Ck+iiCi)\0) (10.1) 

only if the numbers m €z N and k Z are chosen as m = n and k = —1, where n is the index 
of the hamiltonian. In this case ( |10.1| ) takes the form 

= C„_i(e„)C„_2(en-i) . . . Ci(e2)C7o(ei)|0), (10.2) 
and we can use the "quasi" analogues of relations ( ^.9]) 



An{X)Cn-l{Cn) = -^C„_i (A)^„_i (^n) + ^-^^C„_i (C„) A„_i ( A) , 
A - ?n A - 4„ 

Dn{X)Cn^l{Cn) = (A)D„_i (^n) + (e„)Z)„_i (A) (10.3) 

A - ?n A - 4„ 



(which can be easily obtained from basic relations ( |9.1| )) to transfer the operators A„(A) and 
Dn{X) to the right. Exactly as in the case of quasi-Gaudin algebra, each permutation of A- 
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or D-operators with C-operators forming Bethe vector decreases the indices of the former 
by one. Finally, when all permutations will be completed and the operators j4„(A), Z?„(A) 
together with the daughter operators An{^i) , Dn{Ci) , S, = l,...,n appear in front of the 
lowest weight vector |0), they will have index and therefore can be absorbed by this vector 
by means of formulas ( |9.6D . The condition of the cancellation of unwanted terms leads to the 
sysem of Bethe ansatz equations 

n |i:iff^ = 5|', , = i,...,„, (10.4) 



whose form exactly coincides with (^.10), and the corresponding eigenvalues read 



m = a{X) n + d{X) n ^y^, (10.5) 

i=i ^» i=i ^ ^» 

which is again exactly the same expression as in Summarising, we can say, that the 

use of quasi-3^[sZ(2)] algebra enables one to construct an infinite set of completely integrable 
systems, each of which is only partially solvable and has only a certain finite number of 
solutions (of course, provided that both a(A) and h{X) are rational functions). The second 
interesting fact is that the we obtained exactly the same set of solutions as in the case of 
ordinary Yangian. But now these solutions are distributed between different models. 



11 The existence of quasi- Yangian 

In this section we present one of possible realizations of quasi- Yangian. In full accordance 
with the Gaudin case, we construct it from generators of the ordinary Yangian. Assume 
that the latter is characterized by four generators A{X),B{X)^C{\),D{X) acting in a certain 
representation space with lowest weights a?{X) and d?{X). Introduce the limiting operator 

5° = ^ hm X{A{X)-D{X)) (11.1) 



and its lowest weight 



fO = — lim X{a\X) - cf{X)). (11.2) 

2r] A^oo 



Then the needed realization reads 

A^{X) = A{X) ( 1 + r/ 1 + B{X) U ^— 1 , 

B^iX) = B{X) (^1 - , j - A{X) j , 

a(A) ^ C(A) (l + r-''tr/^'^' ) + ^(^) ) ' 

Z.„(A) ^ n^X) (l - ,E1^^1±^±1I1±^) _ c(A) (.^^^) . (11.3) 
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One can easily check by direct calculations that these operators actually satisfy the quasi- 
commutation relations (9J-) and thus form the quasi-3^[sZ(2)]. Since the representation space 
of this quasi- Yangian coincides by construction with the representation space of the ordinary 
Yangian, we can check formulas (|9.6D by using their Yangian analogues (^]^). Taking into 
account that 

50|0)=fO|0), (11.4) 

and using formulas ( |8.6| ), we obtain: 

^.(A)|0) = aO(A)(l + ry!^^i/^)|0), 

D^m = A^) (l - ^ ) |0) - (r^j^) g(A)|0), 

= -aO(A) (7?^) |0). (11.5) 

We see that the lowest weight vector is actually annihilated by "lowering" operators Bn{X) 
and is an eigenvalue of operators Dn{X) only for n = 0. Taking n = in ( |ll.5| ) we obtain the 
expressions for the lowest weights of quasi- Yangian 

a(A) = aO(A)(l + ^^^), <i(A) = ^A) (l - . (11.6) 



Now it is absolutely clear that algebra (|9.lD is a continuous deformation of the Yangian (8.2). 
The role of the deformation parameter is played by o. If a ^ 00, then the n-dependence of 
the operators ( 11. 3| ) disappears and they transform into the ordinary generators of Yangian. 



Respectively, the quasi-relations in (9.1) become the ordinary ones. As to the formulas (|9 



defining the representations of algebra (9.1), they, in the limit a ^ 00 also transform into 
the relations (|8.6|) for the Yangian. 



12 Quasi-exact solvability of quasi-XXX models 

Substituting (pT^ ) into (p^^), we obtain an explicit form of hamiltonians of quasi-XXX 
models: 

^n(A) = ^;"~'^ (A(A)+i^(A)) 
A — a 

+ -^(^(A)-Z)(A))(F0-50 + n + 6) + 
A — a 

+ -^5(A)(F°-S° + n + 26)- 
A — a 

- -^C(A)(F°-5° + n). (12.1) 



Let us now check in an independent way that the operators ( 12.1 ) actually describe quasi- 
exactly solvable models. Denote by the linear hull of all vectors C(^fc) . . . C(^i)|0) with 
arbitrary and k < n. It is known that dim<I>„ < 00 for any n if a^{\) and 

(P{X) are rational functions. From the obvious relations B{X)^n C ^n-i^ C{X)^n C ^n+i, 
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A(A)$„ C D(A)$„ C ^-n, (FO - 5° + n)$„ C and (F^ - 5° + n)$^ C $™ for 

m 7^ n it immediately follows that each of the operators Hn{\) admits only one algebraically 
constructable invariant subspace, This subspace is finite-dimensional and therefore the 



models (12.1) are quasi-exactly solvable. 



13 Historical roots 



Quasi-exactly solvable models came to the light about ten years ago. After first examples 
which have been found heuristically ||, |3^, the three general theories explaining the 
phenomenon of quasi-exact solvability and presenting constructive ways for building such 
models have been formulated. These theories are known under names of "partial algebraiza- 
tion method" |2^, 23, 14, 15, 26, 16 1, "the inverse method of separation of variables" 
|41, 46| and "the projection method" [^, 46]. In sections 14 - ^ we try to interpret the 



results obtained in previous sections from the point of view of these three methods and show 
that, in contrast with the first method, which is not applicable to the study of quasi- Tr alge- 
bras at all, the last two ones enable one to reproduce some general features of the proposed 
formalizm. Especially, this relates to the third, projection method, by means of which we 
essentially obtained all the results exposed in sections |2| - 



14 The partial algebraization method 



One of the first theories of quasi-exactly solvable systems proposed in 1988 in ref. [37| and 
called the "partial algebraization method" was based on a purely algebraic construction. Due 
to simplicity of this theory, its essence can be formulated in several words. 



The method 

The main idea of the partial algebraization method was based on the observation that any 
finite-dimensional representation of any Lie or Lie^ algebra can be realized by means of 
first order differential or pseudo-differential operators acting in the spaces of functions. Let 



Li 



1, 



, d be such generators. Taking their combinations of the form 



-f^ — ^ CikLiLk + ^ CiLi 



;i4.i) 



i,k=l 



1=1 



one obtains a certain second-order differential or pseudo-differential operator having (by 
construction) a finite-dimensional invariant subspace. This means that for this operator at 
least a certain finite part of the spectrum can be constructed algebraically (for more detail 
see refs. 



Note also that this approach is easily extendable to 




the matrix case (see e.g. refs. 



Discussion 

The partial algebraization method seems to be rather general. The models obtainable in such 
a way are not-neccessarily integrable because only an existence of algebraically constructable 
finite-dimensional invariant subspaces is required. But just because of the generality of this 
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method, it contains an essential deficiency not giving a constructive way of solving the re- 
sulting algebraic problems. Essentially, the method stops as soon as the representatibility of 
a certain interesting hamiltonian in the form ( 14. 1| ) is proven. 



What can we say of quasi-Gaudin and quasi-XXX spin models from the point of view of 
this theory? Unfortunately, nothing, because neither the generators of the initial algebra 



used in formulas ( |2.3D and (|8.3D nor their quasi-couterparts used in (|3.3|) and ( |9.3|) are assumed 
to realize any finite-dimensional representation. The problem of finding algebras and their 
finite-dimensional representations which would allow one to represent our hamiltonians in 
the form ( |14.1| ) is interesting and we hope that it will find its solution in the nearest future. 

15 The inverse method of separation of variables 



Another oldest theory of quasi-exactly solvable systems proposed in 1988 in ref. |41] and 
developed in ||4^, |4^ was purely analytic. It can be divided into two parts: 1) construction of 
one-dimensional exactly- and quasi-exactly solvable multi-parameter spectral equations and 
2) reduction of these equations to ordinary one-parameter exactly- or quasi-exactly solvable 
spectral equations in a multi-dimensional space. Consider this method for the simplest case 
of second-order equations. 

The first part of the method 

Let d denote the differential operator, 90(A) = </)'(A), and let F{\) be a given rational 
function. Consider the equation 

{d-F{X))H{\) = E{\)^{X) (15.1) 

for two functions E{\) and ^'(A). It can be shown that it admits a very simple class of 
solutions 

n 

^(A) = n(^-^^)' « = 0,1,2,..., (15.2) 
1=1 

E{X) = F'{X) + F^(A) + 2 X: ~_ f , n = 0,l,2,..., (15.3) 

i=l ^* 

where, for any given n, the numbers Ci, . . . , satisfy the system of numerical equations 

n ^ 

J2 F3^ + ^te) = 0' i = l,---,n. (15.4) 

Before discussing the mathematical meaning of the obtained solutions (which obviously form 
only a very small subset in the set of all possible solutions of equation ( |15.lD ), let us consider 
an analogous construction for pseudo-differential operators. Let denote the shift operator, 
Tfj(j){X) = (j){X + rj), and a(A),d(A) be some rational functions on A. Consider the equation 

(a(A)r^ + d(A)r-i)vI/(A) = E{X)^{X) (15.5) 
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for two functions E{X) and *I'(A). It can be shown that it admits a very simple class of 
solutions 

n 

M/(A) = n(A-ei), n = 0,1,2,..., (15.6) 

1=1 



E{X) = a{X)f[ \^_^^/^ +diX)f[ \^\ \ n = 0,l,2,..., 



where, for any given n, the numbers Cii • • • > satisfy the system of numerical equations 



15.7) 



1, . . . ,n. 



:i5.^ 



An analogous construction for pseudo-differential equations of a little bit different type was 
considered recently in ref. pg|]. 

Now let us explain the reason for which these solutions are interesting to us. A remarkable 
fact (distinguishing these solutions from an infinite number of other solutions) is that for any 
rational function F{X) (in differential case) and for any pair of rational functions a(A) and 
d{X) (in pseudo-differential case) the function E{X) is representable in the form 



N 



E{X) = Wq{X) + E^Wa{X), 



:i5.9) 



where is some number depending only on the form of function F{X) or, respectively, of 
a(A) and d(A), and Wn{X) are fixed functions not depending on ^j. All the dependence of 
E{X) on numbers satisfying the equations ( 15.4 ) or ( 15.^ ) is concentrated in the numerical 
coefficients Ea- It can be shown that for any given n both equatio n (|15.4| ) and ( |15.8| ) have 
exactly {N — 1 -|- n)\/{{N — l)!)n!) solutions. Since n in formulas (|15.2|), ([l5l^) and ( |l5l6D , 
(15.7) can be made arbitrarily large, a total number of algebraic solutions of equations ( |15.1 ) 
and (|15.5| ) is infinite. These reasonings enable one to interpret ( 15.1 ) and ( 15.5 ) as a cer- 
tain exactly (algebraically) solvable A-parameter spectral equation. Indeed, introducing the 
notation 



D{X) 



in differential case, 



{d-F{X)f-Wo{X), 
a{X)T.ri + d{X)T~^ — Wo (A), in pseudo-differential case, 



(15.10) 



we can rewrite the equations ( |15.1| ) and ( |15.5D in the following common form 



N 



DiX)^iX)= [J2EaWaiX)]^{X). 



(15.11) 



\a=l 



In particular case, when A = 1, (15.11) reduces to an ordinary one-parameter spectral 
equation which, in turn, can easily be reduced to the Schrodinger form. This gives us the list 
of known one-dimensional exactly solvable models. 



A remarkable feature of equation (15.11) is that the spectrum of one of its spectral 
parameters is always degenerate, i.e. depends only on n but not on the numbers ^j. Assume 
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for definiteness that the parameter Ejsf has such a degenerate spectrum. Then, introducing 
the notation En = En,N and taking 

D (X) = l ^^^^^^ ~ " ^n,NWN{X), in differential case, 

" I a{X)Tfj + d{X)T,^^ — Wo{X) — En,NW]\f{X), in pseudo-differential case, 

we can represent (|15.11D as an infinite sequence of {N — l)-parameter spectral equations 

L>„(A)^'(A) = ^,VK,(A)j ^(A), n = 0,1,2,... (15.13) 

each of which will obviously be only quasi-exactly solvable. For N = 2 these equations can 
again be considered as one-parameter spectral equations easily reducible to the Schrodinger 
form. This gives us the list of all known series of one-dimensional quasi-exactly solvable 
models. This completes the first part of the method. 

The second part of the method 

Now, what to do if > 1 and > 2 in the first and second cases, respectively? Remember 
that any one-dimensional d-parameter spectral equation with spectral parameters Ei, . . . , Ed 
can be interpreted as a result of separation of variables in a certain d-dimensional separable 
(and hence integrable) system. The parameters Ei, . . . , E^ play in this case the role of 
separation constants and their admissible values have the meaning of the eigenvalues of some 
commuting d-dimensional operators — the integrals of motion of the abovementioned system. 
Note that the form of these integrals can easily be reconstructed by means of the so-called 
inverse procedure of separation of variables. 



Let us first find the N integrals of motion. Hi, . . . , i^^v for the first equation ( 15.11 ) and, 
after this, construct their combination 



N 



H{X) = WoiX) + J2H^Wa{X). (15.14) 



a=l 



The operators H{X) will obviously commute with each other, 

[H{X),H{fi)]=0, (15.15) 



and will have the eigenvalues given by formulas ( |15.3 ) or ( |15.7| ). This means that the system 



described by these operators will not be only integrable but also exactly solvable. 

Let us now find the integrals of motion associated with the second system of multi- 



parameter spectral equations ( |15.13 ). Since the form of these equations depends explicitly 
on n, the obtained integrals also should depend on n. Denoting them by Hn,i, • • • , Hn^N-i, 
let us construct their analogous combinations 



Af-l 



Hn{X) = Wo{X)+En,NWNiX)+Y.Hn,aWaiX), n = 0, 1, 2, . . . (15.16) 



a=l 



From this formula it immediately follows that operators Hn{X) form again commutative 
families 

[Hn{X),HM]=0, (15.17) 
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for any n, but they are not obliged to commute with each other if their indices do not coincide. 
The eigenvalues of operators Hn{\) will be described by the same formulas ( 15. 3| ) or ( |15.7 ) 



but with fixed n. Therefore all the integrable models represented by these operators will 
have only a finite number of explicitly constructible solutions and thus will be quasi-exactly 
solvable. This completes the second part of the method. 

Discussion 

Now what can we say of this method from the point of view of results obtained in sections 



g - |12[ ? A little bit more than in the case of partial algebraization method. First of all, we 
arrived to the same situation as in these sections. The quasi-exactly solvable models have one 
and the same set of solutions as the exactly solvable ones but these solutions are distributed 



between different models. Second, the spectra ( |15.3| ) (resp. (15.7)) of our models and the 



form of auxiliary conditions ( |15.4D (resp. ( |15.8| )) for numbers exactly coincide with the 
spectra and Bethe ansatz equations for the Gaudin (resp. XXX) and quasi-Gaudin (resp. 
q\xas\-XXX) models. This suggests that we obtained the same models but in a very special 
coordinate form. 

Summarizing, one can say that the inverse method of separation of variables correctly 
reproduces some external featurs of models discussed in sections ^ - 12 . However, and 



this is quite obvious, the analytic formalizm exposed in this section is not convenient for 
investigating the properties of the underlying ^[s/(2)] (resp. 3^[s/(2)]) and quasi-^[s/(2)] 
(resp. quasi-3^[s/(2)]) algebras. Moreover, it is applicable only to the study of separable 
models and thus cannot pretend on generality. 



16 The projection method 

The third method of building quasi-exactly solvable models proposed in 1992 in ref. B 



and developed in [46| is the so-called projection method. This method is applicable to all 
integrable and exactly solvable systems with global symmetry. The idea of this method is very 
simple and, in fact, it is implicitly encoded in the inverse method of separation of variables 
which we discussed in the previous section. Indeed, as it follows from the above discussion, 
the complete separation of variables in both A^-dimensional exactly solvable and {N — 1)- 
dimensional quasi-exactly solvable models (constructable in the framework of the inverse 
method) leads (by construction) to one and the same multi-parameter spectral equation. 
But this suggests that the {N — l)-dimensional quasi-exactly solvable models can be obtained 
from their A'^-dimensional exactly solvable conterparts by means of the partial separation of 
variables. Why this observation is so important for us? Because it gives us a possibility of 
formulating a new more general method for constructing quasi-exactly solvable models. The 
generality of this method follows from the fact that it does not require anymore a complete 
separability of variables in the initial exactly solvable model. The only thing which one 
needs indeed is the requirement of a partial separability which, roughly speaking, can always 
be satisfied if the hamiltonian of an initial exactly solvable model admits some global Lie 
symmetry. Let us now explain in more detail what we mean. 
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The method 

Let H{X) be some quantum hamiltonian acting in a certain linear vector space <1>, and let the 
spectral problem 

H{X)(t) = E{\)cp, (16.1) 

have an infinite and discrete set of exactly constructable solutions. Assume that it is possible 
to represent H{\) in the form of a differential or pseudo-differential operator in certain 
variables xi, . . . , and consider $ as a space of polynomials in xi, . . . , x^. Obviously, for 
most of interesting systems this is possible. 

Assume now that there exists some operator L which: a) commutes with H{\): 

[H{X),L\=Q, (16.2) 

b) is representable in the form of a first-order differential operator in the same variables 
Xl, . . . ,XAr: 

N Q 

L = ^ Aj(xi, . . . ■,xn)-q— + B{xi,. . .,xn) (16.3) 

i=l * 

and c) has an infinite and discrete spectrum in <I>. Again, for most of interesting systems 
such operator does actually exist. 

Let In, n = 0,1,2... and n = 0,1,2,... denote the eigenvalues and corresponding 
eigensubspaces of the operator L in <I>. It follows from (16.2) that all these eigensubspaces, 
which, by definition consist of polynomials in xi, . . . ,xn, are simultaneously the invariant 
subspaces for H(X). This means that all spectral problems 

H{X)^ = E{X)(t>, n = 0,1,2,... (16.4) 

are well defined and are exactly solvable. 

Let us now denote by the sets of all (not neccessarily normalizable) solutions of the 
equation 

Lij = ln'ip, n = 0,1,2,.... (16.5) 

Due to the absence of the condition of polynomiality of these solutions, these sets should be 
wider than 

^nC^n, n = 0,1,2,..., (16.6) 

but, as before, they will be the invariant subspaces for H{X), and threfore the extended 
versions of equations (16.4) 

H{X)'4) = E{X)i^, ip(^^n, n = 0,1,2,... (16.7) 

will also be well defined and, formally, will be quasi-exactly solvable. This, however, will not 
be that kind of quasi-exact solvability which we discussed before and which we want to have. 
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This will be rather some pathologically trivial case of it, because the hamiltonian H{X) still 
does know that it should act only in one of the sets and acts equally well in all of them. 
Practically this means that we cannot interpret the equations ( p.6.7 ) as spectral equations for 
different models. 

In order to improve the situation, one should restrict somehow the action of the hamil- 
tonian H{X) to one of the sets ^n- But for this it is neccesary to have a little bit more 
information of the structure of these sets. We can extract it from the equation ( [16.5 ) which. 



according to ( 16.3 ), can be rewritten as 
N 



\^^Ai{xi, . . ■^^n)-^J ipixi, . . .,xn) = (^In - B{xi, . . . ,XN)^ipixi, . . .,Xn)- (16.8) 

The general solution of this equation has the form 

ip{xi, . . .,xn) = Pn{xi, . . . ,Xn)f[cri{xi, . . . ,x„), . . .,aN-i{xi,. . .,xn)] (16.9) 



where pn is a certain arbitrarily fixed particular solution of ( 16.81 ) and cjj, i = 1, . . . , N — 1 



are some functionally independent solutions of the auxiliary equation 



/ ^ d \ 

Ai{xi, ^N)-Q^j (^{xi, . . . , xn) = 0. (16.10) 



The function f{ai, . . . , cr^^i) in ( 16. S| ) is arbitrary and therefore all the information of the spe- 



cific properties of the spaces is concentrated in p„. The set of all functions /(<ti, . . . , cjAr_i), 
which we denote by F, contains for any given n a certain subset Fn corresponding to explicit 
solutions of the initial equations (16.4). 



Let us now note that, due to the functional independence of functions 
P = [pn{xi,...,Xn)]", 

ai = ai{xi,...,Xn), i = 1, . . . ,iV - 1, (16.11) 

they, for any given n, can be considered as the new independent variables (instead of 
xi, . . . ,xn). Rewriting the hamiltonian H{\) in terms of these new variables and using 
the obtained expressions ( |lO| ) for funct ions ij) with given n, we can rewrite equations ( |16.7D 
in the form: 

^„(A)/(ai,...,ajv-i) = ^(A)/(ai,...,aiv-i), / e F, n = 0,1,2,... (16.12) 

where 

Hn{\) = p-^H{\)pn, n = 0,1,2,... (16.13) 



Now the only one last step remains. The step which just suggested us to call this method 
the projection method. Let us look at formula ( |16.12| ). We see that the transformed hamil- 
tonians Hn{\) are, by construction, certain differential or pseudo-differential operator in all 
variables p and cxi, . . . , a^-i- At the same time, they act only on functions on A^ — 1 variables 
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cJi, . . . ,o"Ar_i and give again the functions of the same variables. This means that the most 
general form of these operators is 



^„(A)=F„(A) + K(A)^ 



n 



0,1,2,..., 



(16.14) 



where H'^{X) are certain operators in all variables pn and cJi, . . . , a^, and Hn{X) are operators 
in variables ai, . . . , ajy only. This enables one to rewrite (16.12) in the following final form 



Hn{X)fi(Ji,...,aN-i) =S(A)/(fTi,...,cJ7v-i), / e F, 



n 



0,1,2, 



(16.15) 



We see that these equations have the form of ordinary quasi-exactly solvable ones, because the 
set F contains by construction the subsets -F„ in which these equations are exactly solvable. 
All the information of a quasi-exactly solvable model is encoded now in its hamiltonian which 
explicitly depends on n. 



Discussion 

Let us now try to explain the results of sections |2| — 12 from the point of view of the 
projection method. For this, let us assume that H{X) is a hamiltonian of Gaudin or XXX 
model expressed respectively via generators of Q[sl{2)] or 3^[s/(2)] algebras by formulas ( |2.3D 
or (|8.3D . It is known that generators of both C/[s/(2)] and y[sl{2)] algebras can be constructed 
from generators S~, S^, S^, a = 1, . . . , N of algebra sl{2) (g) . . . (g) sl{2) (N times). The 
corresponding formulas have the form 



50(A) 
+S+{X) 



for Gaudin algebra C/[s/(2)], and 

( AW B{X) 
[ C(A) DiX) 




N 

E 



qO 




X — aa 



(16.16) 





(16.17) 



for the Yangian 3^[sZ(2)]. Using formulas ( |2.3| ) and (|8.3D , we can respectively express the 
hamiltonians H{X) of both Gaudin and XXX models via generators S~, S^, S^, a = 
1,. . . ,N. Remember that these generators can be realized as differential operators 



d 

dXa 



_d_ 

' dXr 



_d_ 



(16.18) 



acting on polynomials in variables xi, 



,XN- 



Formulas ( 16.16| ) and ( |16.17 ) together with 



( 16.18 ) enable one to represent the hamiltonian H{X) in the form of some differential operator 
in N variables xi, . . . ,X]\f. The space $ in which this hamiltonian acts becomes respectively 
the space of polynomials. 

Now note that, the hamiltonian H{X) has a global sl{2) symmetry, which, in particular, 
implies its commutativity with the operator of the z-projection of total spin 



L 



N 



(16.19) 



a=l 
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According to ( 16.18 ), the differential form of this operator is 

N Q 

^=E^"5^ + ^' (16.20) 

a=l °' 

where F = Fi + . . . + Fjy. The eigenvalues of this operator in $ are given by = n + F and 
the corresponding eigensubspaces are formed by all homogeneous polynomials of degree 
N. 

To construct the spaces , one should find a general solution of the equation 
/ ^ d \ 

yl2 ^"^^ + ^ ) ^(^^' ■ ■ ■ ' = + • • • , xn) (16.21) 

which, obviously, reads 

. . . , :r^) = x^^f (^,..., ^) . (16.22) 

\XN XN J 

The new variables can hence be chosen in the form: 

P = xn, ai = ^^, i = l,...,N — l. (16.23) 

XN 

According to general prescriptions of the projection method, the first (preliminary) step in 
the reduction procedure lies in rewriting the hamiltonian H(X) in new variables ( |16.23D and 
in constructing the homogeneously transformed hamiltonians Hn{X), according to formula 



(16.13). In our special case this formula gives: 

HniX) = p-''H(X)p^. (16.24) 



Using formulas (|2.3|) and (|8.3|), we obtain for Hn{X): 

H^{\) = Sl{\)Sl{\) - \s-{\)~St{\) - ^5+(A)S-(A), (16.25) 

5°(A) = p-^S\X)p^, S-{\) = p-"5-(A)p-, 5+(A) = p-^S+{X)p^, (16.26) 
for the Gaudin case, and 

Hn{\) = An{\) + Dn{X), (16.27) 

i„(A) = p-M(A)p", 5„(A) = p-"I)(A)p", (16.28) 

for the case of Yangian. An explicit computation of operators staying in the right hand sides 
of formulas (|16.25D and (|16.27[) gives 



/ 50(A) -5-(A)\ / S\X) -p-'S-{X) 

[ +5+(A) -50(A) J [ +pS+{X) -50(A) 

P^+FN+n-a -1;-^ \ 1 
,2 d 



+ 



p^d_ + 2FNP+{n-a)p -p§-^ - Fn - {n - a) ) \ - aN^ ^ ' ' 
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for Gaudin algebra C/[s/(2)], and 



i„(A) BniX) 
CniX) Dn{\) 



A{\) p-^B{\) 
pC{\) D{\) 



1 + 



P§-n + FN + {n- a) 



X-aN \ P^-§-p + 2Fnp + (n - a)p 



d n—c 

"5p p 



'd-p 



Fn - in- a) 



for the Yangian 3^[sZ(2)]. Here 



iV-l 



(J =^aa 



a=l 



d 
dan 



(16.30) 



(16.31) 



and the entries of the matrices stayng in the right hand sides of ( 16.29 ) and ( 16.30 ) are given 
by 



and 



S^iX) -S'{X) 
+5+(A) -50(A) 

/ ^(A) B{X) \ ^ 
[ C{X) D{X) i 



N-l 

E 



rr -I- F --^ 



1 



„2 d _l_ nip rr —rr 

Q=l V ^ada^ -r^ra(Ja "ada^ 



X- an 



N-l 

n 

a=l 



1 + 



V 



rr _|_ 7? d_ 



X — an \ a 



2 d 
a dcTa 



+ 2Fn(Jn 



a 



Fn 



(16.32) 



(16.33) 



In order to complete the reduction of the initial hamiltonian H(X) to a sequence of 
hamiltonians Hn{X) of quasi-exactly solvable models, one should project the operators ( 16.25 ) 
and ( 16.27 ) onto a space of functions depending on variables di, . . . , a^-i only. After some 
simple algebra one obtains 



Hn{X) = 50(A)50(A) - ^S-_,,{X)S+{X) - ^S+^,{X)S-{X) 



(16.34) 



with 



^n(A) 
5- (A) 



S'{X) + 
S%X) + 
S-{X) + 



n — a 
X — Cat ' 
n — a + Fjy 

A — UN ' 
n-a + 2Fn 



for the Gaudin case, and 



A — Oat 
Hn{X) = An{X) + Dn{X) 



(16.35) 
(16.36) 



with 



Bn{X) = B{X)[l-'n 

Cn{X)=C{X) fl + 77 



^(A) 7? 



A — ctv 
n — a 



AniX) = A{X) ( 1 + r? j + B{X) ( r] ^ -— ) 

n- a + 1/2 + Fat 

A — Cat 

n-fj- 1/2 + Fat 



Dn{X) = D{X) 1-r? 



A — aj^ 

n - a + 1/2 + Fn 



+ Z)(A) 77 



A — CAf 

n- a + 2Fn \ 



X — On 



C{X) r? 



A — qn 
n — a 



X — ajsf 



(16.37) 
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for the case of Yangian. 

It is not difficult to see that formulas ( |16.3^ ) and ( 16.36| ) exactly coincide with expressions 



0|) a nd ( |9.3|) for hamiltonians of quasi-Gaudin and quasi-XXX models. As to the formulas 
16.351) and ( |16.37| ), they, up to a trivial change of notations, are identical to formulas (|57 



and ( 11.3 ) for generators of G[sl{2)] and y[sl{2)] algebras. 

Summarizing, we can say that the projection method enables one to reproduce correctly 
the structure of quasi-Gaudin and quasi-XXX models together with all properties of the 
undelying Q[sl{2)] and J^[sZ(2)] algebras. However, this method contains so many technicali- 
ties and is so much model dependent that its practical usefulness seems rather doubtful. In 
next sections I will expose another simple method which practically without any calculations 
allows one to construct quasi-analogue for any algebra. 



17 Graded associative algebras 

Now we start the exposition of a general and model independent method for building quasi- 
Tr algebras and their "representations" . The idea of this method is based on the observation 
that both Tfj and quasi- Tr algebras can be considered as two different realizations of one 
and the same abstract Z^'-gradedpl unital associative algebra A. The difference between 
the "ordinary" and "quasi" versions of this algebra lies in a different nature of its elements 
and also in different definition of the corresponding associative product. In this section we 
remind the reader some general properties of an abstract algebra A realized in an abstract 
linear vector space V . The structure of the pair A, V is described by the following standard 
axioms. 

Axiom 1. The set ^ is a union of non- intersecting layers A^, n € Z'^ each of which forms 
a linear space over C. If a G A"^ then we say that a has grading n and write this as g{a) = n. 
The summation of elements belonging to different layers is forbidden. 

Axiom 2. For any n,m (z Z''' there exists a multiplication map (A",^'") ^n+m ^^^^[q)^ 
puts into correspondence to any pair of elements a G A"' and b G A"^ their product ab G A^~^"^. 
This product is a) associative: a{bc) = {ab)c, for any a G A^, b G A"^, c £ A'- and n,m,l G Z^, 

b) distributive: {a + b)c = ac + bc, c{a-\-b) = ca + cb for any a, 6 G A"-, c G A^ and n,m G Z^, 

c) has a unit I G A^ such that la = al = a for any a G A"' and n G Z^ . 

Axiom 3. The set y is a union of non- intersecting layers V"', n G Z^ each of which forms 
a linear space over C. If u G V"" then we say that v has grading n and write this as g{v) = n. 
The summation of elements belonging to different layers is forbidden. 

Axiom 4. For any n, m G Z"^ there exists a multiplication map {A^,V^) yn+m 
which puts into correspondence to any pair of elements a G A"^ and v G their product 
av G 1/"+™. This product is a) associative: a{bv) = {ab)v, for any a G A"', b G A^, v G 
and n,m,l G Z^, b) distributive: (a -|- b)v = av + bv, a{v + u) = av + au for any a,b G A^, 
v,u G and n,m G Z"^ , c) has the same unit I G Aq such that Iv = vl = v for any v G V"' 
and n G Z^. 

Any A-algebra contains two important subalgebras. One, which we denote by ^4^^, is 
formed by all zero-graded elements of A. The second one which we denote by A~^ consists 

^^Here — Z ® . . . ^ Z {r times) denotes the set of all r- dimensional vectors with integer components. 
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of all elements with non-negative gradings. The non- negativity of a grading means that all 
components of the grading vector n = (ni, . . . , n^) are non- negative integers. The complement 
of A'^ in A, consisting of elements whose gradings contain at least one negative component 
(we call such gradings quasi- negative) , does not form any subalgebra of A. Nevertheless, this 
complement plays an important role in applications and we shall denote it by A~ . 

The space V also contains an important subspace formed by all vectors with nonnegative 
gradings. We denote it by V~^. The importance of the space lies in the fact that it, under 
some special conditions, takes the meaning of an algebraically constructable representation 
space of algebra A. In order to realize these conditions one needs a special so-called lowest 
weight vector |0) with grading zero which is annihilated by any operator from A~: 

a"IO) = 0, for any a~ e A' (17.1) 

and is an eigenvector of all operators from A^: 

a»=e(a°)|0), for any a° S (17.2) 

Here e{a^) denotes some linear functional of which is usually called the lowest weight. 
Once a vector |0) is constructed, one can construct the space in a purely algebraic way, 
simply by acting on it by various elements of A'^: 

V+= \J a+lO). (17.3) 

a+eA+ 



Note that the operator elements of ^-algebras can sometimes be treated as quantum 
observables and the vectors of the space — as physical states. In this case the grading 
of vectors v from can be considered as a collection of "occupation numbers" for certain 
elementary excitations and the action of operators from A~^ and A~ on these vectors — as a 
process of creation and annihilation of excitations. The lowest weight vector takes respectively 
the meaning of the vaccuum state. If the number of excitations is a conserved quantity, then 
the hamiltonian of a system should neccessarily be a zero-graded element of subalgebra A^. 
Fixing the hamiltonian ^ A^, the vaccuum state \0) eV^ and the space we essentially 
fix a physical system. 

Once a physical system is fixed, one can state the problem of finding all its physical states. 
This leads us to an eigenvalue problem for hamiltonian h^: 

/iV = €(/), (/> G V+. (17.4) 

Note that the knowledge of only general properties of algebra A listed above is, obviously, 
not sufficient to present an explicit solution of problem ( [17. 4 ). In order to solve it, we need an 



information of more specific properties of this algebra which strongly depend on a concrete 
form of its elements and therefore are not encoded in axioms 1-4. In fact, what we actually 
need, is the existence of special bilinear relations between the elements of algebra A. 

Indeed, assume that for given € A^ and each m > there exists an element € 
with grading m for which the relations 
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hold. Here no special requirements to elements € A^, a~. € A" and a^. G A'^ are 
assumed. Acting by this bilinear operator relation on |0) and using formulas ( 17.1| ) and 
( 17. 2| ), we obtain 

= EmfAm, em = e(6^), (/)^ = a+|0), meZ (17.6) 

which gives us infinitely many purely algebraic solutions of the spectral problem ( 117.4] ). 
Vectors (pm = o-ml^) solving this problem are typical Bethe vectors. Despite the seeming 
simplicity of this solution, its actual construction is sometimes very complicated. The most 
non-trivial thing is to find appropriate Bethe vectors for a given hamiltonian. Fortunately, 
for many physically interesting systems this problem has been successfully solved (see e.g. 
sections § and P) . 

18 Modification of an associative product 

Up to now the consideration of A-algebras and their representations was rather general: the 
elements of operator algebra A were considered as elementary objects not having any internal 
structure. In this section we consider a special realization of algebra A whose elements are 
some composite operators and their product is defined in a non-standard way. We demon- 
strate that such a realization also has a lot of interesting physical applications. In order to 
distinguish this realization from the general operator algebra A discussed in section ^ we 
shall use for it the notation A and reserve the star-symbol for the product of its elements. 

We define the internal structure of the elements d E A as follows. We consider them as 
functions on Z^ whose particular values a[n] are some linear operators on V. The product of 
two elements d and b will be denoted by d * 5 and defined as 

(d * b) [n] =a[n + g{b)] ■ b[n] (18.1) 

where the dot means the standard operator multiplication. The action of an element d on 
the vector v &V is defined as 

av = a[g{v)]v. (18.2) 

It is not difficult to check that these definitions do not contradict the axioms 1-4. 

As in general case, we can introduce the subsets A^, A'^ and A^ of algebra A consisting 
of the elements with zero, non-negative and quasi-negative gradings. The values of these 
subsets on Z"^ we respectively denote by ^'^[n], A"'"[n] and j4~[n]. 

In order to construct the lowest weight representation of algebra A we need the lowest 
weight vector |0) G y satisfying the hat- analogues of equations ( |17.1| ) and (17.2): 



a-|0) = 0, for any d~ G i", (18.3) 

a>) = e(d°)|0), for any d° € i°. (18.4) 
If such a vector does exist, we can define the representation space as 

V+= [j d+|0). (18.5) 
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Let us now rewrite formulas ( 18.3 ), ( 18.4 ) and ( 18.5 ) in components. Using definition ( |18.2 ) 
and taking into account tliat the grading of the lowest weight vector is zero, we obtain: 



a-[0]|0)=0, for any a-[0] G ^"[0], (18.6) 

a°[0]|0) = e(a°[0])|0), for any a°[0] G A°[0], (18.7) 

and 

V+ = y a+[0]IO). (18.8) 

a+[0]eA+[0] 

From the above consideration it follows that the lowest weight vector |0) is not obliged to 
be annihilated by all "lowering" operators a~ [n] (with n 7^ 0) and to be an eigenvector of all 
"neutral" operators aP[n] (with n ^ 0). 

Why the yl-algebras maight be interesting to physicists? There are several reasons for 
this. Let us try to explain one of them. Assume that we managed to realize one and the same 
abstract ^''-graded unital associative algebra (abstract in the sense that it is specified only by 
some relations between its elements without any concretization of their internal structure) in 
two different ways: 1) as an operator algebra A and 2) as an algebra of composite operators 
A, i.e. in terms of the dot- and star-products, respectively. Then we can say that algebras A 
and A are isomorphic. 

Assume now that algebra A contains two operators a and b commuting with each other: 

a-h = b-a. (18.9) 

If we consider this relation from the physical point of view, we can imagine that the operators 
a and h represent two observables which can be measured simultaneously. After this one can 
try to make some conclusions about their common spectral properties, etc. In other words, 
the relation ( 18. 9| ) may be quite meaningful and informative for physicists. 



Due to the assumed isomorphism between the A and A algebras, an analogous relation 
can also be written for the corresponding composite operators a and h. They also should 
commute with each other: 

a*6 = 6*a. (18.10) 



But what can we now say of formula ( 18.10| ) as physicists? At first sight, nothing, because 



the composite operators and their star products do not have themselves any clear physical 
meaning. Remember, however, that the composite operators have the meaning of the sets of 
ordinary operators and their star product can be expressed via ordinary operator product. 
May be, the relation ( |18.1[1| ) will become more meaningful after rewriting it in components? 
Let us see. Using formula (|18.lD , we can write: 

a[n + g{b)] ■ h[n] = b[n + g{a)] ■ a[n] (18.11) 

Now the relation ( |18.11| ) contains everything what we need: the operators and operator prod- 
ucts. However, we also see that this formula does not describe any longer a commutativity 
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of any two objects: in general, this is a certain relation between four different objects! Let 
us assume, however, that both o and b elements have zero grading. In this case the situation 
drastically changes: the relation ( |18.11| ) takes the form 

a[n] ■ b[n] = b[n] ■ a[n] (18.12) 

and says us that the operators a[n] and b[n] commute in the standard operator sense for any n. 
So, we see that we obtained even more than expected. Instead of a single pair of commuting 
operators which we started with, we obtained an infinite set of pairs of such operators! 

What does this mean from the physical point of view? There are many examples of A- 
algebras containing large families of mutually commuting operators. If the grading of these 
operators is zero, then it is possible to interpret them as integrals of motion (hamiltonians) of 
some completely integrable quantum systems. Assume now that we managed to construct a 
A-analog of the above A-algebra which is isomorphic to A. Then we can repeat the reasonings 
given above and, instead of a single family of commuting operators representing a single 
integrable system, construct an infinite and discrete set of families of commuting operators 
which can be treated as integrals of motion of infinitely many different integrable quantum 
systems. 

If we actually want to treat the infinite sets of commuting operators appearing in our 
scheme as hamiltonians of some physical systems, we should discuss the possible ways of 
solving the corresponding spectral problems. Let us try to do this in the same way as in 
the case of general ^-algebras (see section 0). Let be a zero-graded element of algebra 
A considered as a hamiltonian of some integrable quantum system. Assume that spectral 
problem for this hamiltonian is exactly solvable. This implies the existence of the relations 
( 17. 5| ) in algebra A. Assume also that we have in our disposal an algebra A isomorphic to A. 



Then, using this isomorphism, we can rewrite relations (17.5) in the A-algebraic form: 



^° *a+ = a+ + ^a+^ *a^^. (18.13) 

i 

Here, as before, m is an arbitrary non-negative element of Z*" and no special requirements 
to the elements € a~. € A~ and a+, € are assumed. Acting by ( |18.13 ) on the 
lowest weight vector |0) and using formulas ( |l8^ ) and ( |l8lD we obtain 

/iV„ = e{b^J<Pm, <Pm = a+ |0). (18.14) 



Formula ( 18.14 ) does not have yet any clear physical meaning. In order to make it physicaly 
meaningful, one should rewrite it in components. Taking into account that the grading of 
vector is m and using formula ( 18. 2| ), we obtain immediately 



h"[m]<P^ = em<Pm, em = e(6^[0]), = a+ [0]|0), m e (18.15) 

which gives us a purely algebraic solution of the spectral problems for h^[ni\, m ^ Z"' . Now 
however, this solution is incomplete, because the vectors (pm = (i-tnl^'W^) with any given m 
form only a negligibly small part of the whole space . This enables one to make a general 
conclusion that if the mother integrable system associated with the A-algebra is exactly 
solvable, i.e. admits a complete algebraic solution of the spectral problem, then the daugther 
integrable systems associated with the ^-algebra will be only quasi-exactly solvable, i.e. only 
some finite parts of their spectra can be constructed algebraically. 
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19 Construction of modified graded associative algebras 



In this section we describe a simple scheme which essentially opens a practical way for building 
the graded algebras of composite operators. We demonstrate that any ^-algebra whose 
elements admit a special semi-differential realization can be reduced to a certain ^-algebra 
which is isomorphic to A. 

Assume that each element a of a given algebra A (which is obviously an operator) admits 
the representation 

a = e'-a^''\h{a) + h,{a)di + hikia)d^dk + . . .) (19.1) 

where t = {ti, . . . ,tr} is a certain variable vector and di = d/dti, i = 1, . . . ,r are partial 
differential operators with respect to the components of t. The expression in the exponent is 
a scalar product of two r-dimensional vectors t and g{a). The symbols h(a), hi(a), hn^^a), . . . 
denote some i- independent linear operators acting in a certain vector space W. A concrete 
form of these operators depends on the element a. The summation over repeated indices 
is also assumed. In this case, the elements of A can be viewed as operators in the space 
W = V X T where T denotes the space of all analytic functions of t. Assume also that any 
vector w of the space W admits the representation 

w = e^-9('"K{w) (19.2) 

where v{w) denotes some t-independent vector of V. 

Let us now associate with elements a defined by formula ( |19.1 ) some operator- valued 



functions a on Z*" whose particular values, a[n], are operators in the space V defined by the 
following relations 

a[n] = h{a) + hi{a)ni + hik{a)nink + ... (19.3) 

It is not difficult to check by means of direct computations that the elements a form an A- 
algebra with the products defined by formulas ( |18.1| ) and ( |18.2| ). This means that any algebra 



A whose elements admit the representation ( |19.1| ) and act in the space V whose vectors are 
representable form (19.2) can be reduced to a certain A algebra. 



The next important observation which can be immediately made after analysing formulas 
( 19.1| ),( p^9.2| ) and ( |19.3D is that the map A ^ Awe constructed is an isomorphism. This means 



that if there are some special (polynomial) relations between the elements of the operator 
algebra A then the map A ^ A will induce the same relations between the corresponding 
elements of A. The only difference will concern the definition of the product because the 
ordinary operator product in A should be replaced by a rather specific product in A. 



20 Quasi-Lie and quasi-Lie^ algebras 

It is known that any Lie or Lie^ algebra C of the rank r can be considered as a partially 
Z''-graded algebra. The grading can be introduced in many different ways, but there is a 
special, the so-called root grading, which seems to be most natural. In this case the gradings 
of separate elements of a given algebra form a certain finite subset in which describes the 
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root system of an algebra in the basis of r simple roots. This fact enables one to endow with 
grading also the associated universal enveloping algebra U{C). This, obviously, will be a full 



^''-grading which makes it possible to consider IA{C) as an A-algebra described in section 17 



Now note that each algebra C can be realized in the form of differential or pseudo- 
differential operators in several variables xi, . . . , xat. A minimal possible number N of these 
variables is r (i.e. is equal to the rank of algebra C and to dimension of grading vectors), 
and the maximal one is not limited. The differential operators representing the elements 
of algebra C may depend additionally on certain numerical parameters characterizing the 
representation of algebra C and determining the eigenvalues of zero-graded elements of C 
Therefore, the grading of these parameters also should be equal to zero. From this obser- 
vation it automatically follows that the variables xi, . . . ,xn should neccessarily be graded. 
More exactly, at least r of these variables should be graded and these gradings should be lin- 
early independent. Otherwise, it would be impossible to construct from them the generators 
associated with simple roots whose gradings are linearly independent by definition. 

Taking into account the last reasonings, let us now try to reduce the generators L"" of 
algebra £ to a special canonical form by means of an appropriate change of variables. Note 
that the existence of r variables with linearly independent gradings makes it possible to 
choose the new ones in such a way that to have r variables pi, p2, ■ ■ ■ , Pr with gradings 
{1,0,..., 0}, {0, 1, ... , 0}, . . . , {0, 0, . . . , 1} and all the remaining variables with grading zero. 
The most general form of generators L"" in the new variables reads 

= pf ■ ■ ■ pf [ii^ + Y.Htp~ + ..)i (20.1) 

where = (nf,...,n°) is the gradings of and ff", Hf, etc. are certain differential 
operators in the remaining variables. The expression ( |20.1| ) is, however, rather cumbersome 
and therefore, it is convenient to make one extra change of variables, namely, /7j = e*% i = 
1, . . . , r, which brings ( |20.1 ) to the desired canonical form 



L'' = e^9{L-\H'' + H^di + ...) (20.2) 

with t = {ti, . . . ,tr}- From this expression it immediately follows that each element of the 
universal enveloping algebra U{C) is also representable in the form ( |19.1 ) and for this reason 
U{C) can be reduced to a ^-algebra isomorphic to A. 

Let us now look at the explicit formulas. Let £ be a Lie algebra. Then the commutation 
relations for its generators read 

L" ■ - ■ L" = Cf L\ (20.3) 

Rewritting the star-analogues of these relations (which read similarly) in components L"[n] 
we obtain the relations 

^{n + g{L^)\ ■ L\n\ - L\n + giL")] ■ L°[n] = Cf L^[n] (20.4) 

which do not form any Lie algebra with respect to the ordinary operator multiplication. We 
call such an algebra the quasi-Lie algebra and the expression staying in the left-hand side of 
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(|20l ) — the quasi-commutator. Let now £ be a Lie^ algebra. Then the most general form 
of the relations for its generators is 

L^-L^ = R'J^ • L-^. (20.5) 

The corresponding star- analogues of these relations being rewritten in components will read 

L'^ln + g{L')] ■ L\n] = Rf^ L'[n + g{L^)] ■ L''[n]. (20.6) 

These relations also do not form any Lie^ algebra. We call an algebra described by them the 
quasi-Licg algebra. 

Consider an example. In section ^ we already considered the case of quasi-s/(2). Consider 
now its Q'-deformed version. Let £. = slq{2) with four generators a,b,c,d with gradings 
0,-1, +1,0, respectively. Then the standard defining relations for this algebra read 

ab = qba, 

cd = qdc, 

ac = qca, 

hd = 



with an additional requirement 



ad — da = qbc — q ^cb, 

be = cb, (20.7) 



det 1^ ^ ^^j=ad-qbc = 1. (20.8) 

Using the simplified notations dfi — (i[Ti\^bfi — 6[?T/],C72 — c[?T/],(i72 — (^[tz] , we caiii write down 
the quasi-counterparts of relations (20.7) and ( |20.^ ) which respectively will read 



a-n-lbn 


= qbuCln, 






C-ndn 


= qdn+lCn, 






bndn 
^ndn ^n^n 
bn+lCn 


— qC-nO"ai 
= qdn-lbn, 

= qbn+lCn - ^ 

— C-n—lbn 


Cn~lbni 


(20.9) 


a^dn 


- qbn+lCn = 1- 




(20.10) 



and 



We call an algebra defined by these relations the "quasi-s/g(2)" . 
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21 The main theorem 



Collecting the results of sections 17 - 20 we can formulate the following theorem. 

Theorem. Let Che a certain Lie or licg algebra defined by commutation relations ( ^0.31) 
or ( pp. 5 ). Assume that its generators obey some polynomial relations 

j,k a.-i,...,ak 



n 



derivable from the defining relations ( |20.3| ) or ( 20. 5| ) of algebra C. Then the generators L 
of an abstract quasi->C algebra defined by quasi-commutation relations (20.4) or ( pO.6 ) obey 
the "quasi" analogues of relations (21.1) which read 



j,k aiviifc 



and in which 



50 = 0, gk = Y,g{L^^] 



(21.3) 



1=1 



Proof. Denote by A the universal enveloping algebra of algebra C. It obviously can be 
considered as a Z'^'-graded associative algebra with unit. Using generators L"'[n\ of a given 
quasi->C algebra, we can construct the composite operators L"" whose product is defined by 
formula ( 18. 1| ). The algebra C of these operators is isomorphic to C by construction. Let us 
now construct the universal enveloping algebra A for C. The isomorphism of algebras C and 
C implies the isomorphism of their universal enveloping algebras A and A. This, in turn. 



implies the existence of hat-versions of relations ( 21.1 ) in A. Rewriting these hat-relations in 

which completes the proof. 



the components we obtain formulas (21.2) and ( |21.c 

Corollary 1. All relations in algebra Tr imply the existence of analogous quasi-relations 
in the corresponding quasi- Tr algebra. These quasi-relations can be obtained in the following 
simple way. One should consider separately each monomial in an initial T/j-algebraic relation 
and endow all its co- factors by one and the same index n € Z"' . After this one should shift 
the index of each cofactor by a total grading of all right-standing co- factors. Collecting 
all monomials transformed in such a way together, we obtain the desired form of quasi- T- 
algebraic relation. 

Corollary 2. Each integrable and exactly solvable model associated with algebra Tr 
has an integrable and quasi-exactly solvable conterpart associated with the corresponding 
quasi- Tr algebra. Applying the rules given in corolary 1 to formulas describing the form 
of commuting integrals of motion and determining the structure of Bethe ansatz solutions 
of integrable and exactly solvable models, we can obtain the corresponding formulas for 
integrable and quasi-exactly solvable models. In particular, all formulas of sections 
P - 12 can be obtained elementary from the analogous formulas of sections § and 




From these two corollaries it immediately follows that the only non-triviality in con- 
structing the quasi-exactly solvable modifications of integrable exactly solvable models lies 
in finding the realizations of quasi- Tr algebras defined by quasi-commutation relations 



Rai3-^s{>^ - lJL)T^p{\n + gs^)Tsa{l^,n) = T/3sin,n + ga'r)Ta'yi\n)R^Spai>^ - /^) (21.4) 
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with ga/3 = g{Ta/3{X,n)). One of the possible ways of constructing such reaUzations was 
described in section |l^. It is quite obvious however that the proposed reaUzations do not 
expire the set of all possible ones. This can be seen even from the defining relations (3.1) 



and ( |9.lD for quasi-t/[s/(2)] and quasi-3^[s/g(2)] algebras which leave a considerable freedom 
in choosing the form of their generators (see formulas ( |3.2| )and ( |9.2[ )). The problem of listing 
all non-equivalent realizations of these algebras is only a part of the following general one: 

Problem. Construct the representation theory for quasi- Tr algebras. 

We expect that solution of this problem may lead to revealing new interesting math- 
ematical structures and notions which neccessarily find their applications in the theory of 
quasi-exactly solvable systems. 

22 The classical limit of quasi- Yang-Baxter algebras 

It is known that the Gaudin algebra t/[s/(2)] has a well defined classical conterpart — the 
so-called classical Gaudin algebra. The corresponding Lie-Poisson commutation relations for 
its generators (which in this case take the meaning of functions on the phase space) read 

{50(A), 5° (/x)} = 0, 

{S+{X),S+ip)} = 0, 

{5-(A),5-(/i)} = 0, 

{S\X),S^{^,)} = -(A-/.)-i{5+(A)-5+(/i)}, 

{50(A), 5- (/x)} = +(A-/.)-i{5-(A)-5-(/i)}, 

{5-(A),5+(/x)} = -2(A-/x)-H50(A)-50(/x)} (22.1) 

and imply the commutativity of functions 

H{X) = S°{X)S°{X) - 5+(A)5-(A) (22.2) 

playing the role of integrals of motion of a classical completely integrable system. The latter 
is usually reffered to as the classical Gaudin model. 

It is naturally to ask ourselves if there exists any classical limit of the quasi-Gaudin algebra 
defined in section ^. The answer to this question is positive. In order to demonstrate this 



we note that the quasi-Gaudin algebra defined by quasi-commutation relations (3.1) admits 
a class of equivalent representations which can be obtained by an appropriate rescaling the 
integer parameter n. Indeed, it is not difficult to see that if one uses instead of (|3.1| ) the 
following modified relations: 

Snh-hWSnhif^) - S~f,^f,{^l)S~f^{X) = 0, 

^nn+n(A)5+(//)-5+(/i)50^(A) = -^(A-/.)-H5+(A)-5+(/i)}, 
SL-nWS-nif^) - S-^if,)SU>^) = +h{X - fi)-' {S'^iX) - S-^{f,)} , 
Snn+nWSM - Sin-nif^)S-nW = -2?i(A - fJ^r'{SU>^) - SU^^)}, (22.3) 
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in which ?i is a parameter, then the operators 

HnniX) = 5°,(A)S°,(A) - \s-n+nWS^nW - ^5+ _,(A)5-,(A) (22.4) 

will as before commute with each other for any given n. Let us now rewrite formula ( p2.3D 
in the form 



0, 



[5+(A),5+(m)] , g+;,+JA) - S+ (A) ^^^+.(^1-5^ _ 

h + n ^^^^^^ h ^-^^^^ - °' 

KnW,SM] , JA) - S°,(A) g+(A)-5+(/x) _ 

[SLW,S-n{^^)] , g°,-a(A) - 5°,(A) g-,(A)-g-,(/x) _ 

['S'na(A)) -S'itsC/^)] , 'S'n;i+?i(A) - 'S'„;i(A) + 'S';[7j_;j(^) - S';[;j(^) _ 
^ ^ ^ ^nnyl^) ^ ^nny.^) - 

A — /i 



(22.5) 



and apply to it a standard dequantization procedure. For this one should assume that the 
parameter fi is small, n is large, and their product t = nh is finite. Considering t as a 
new continuous variable and taking the limits ?i — s- and n — > oo, we can replace the 
commutatators in ( ^2.5 ) by Poisson brackets and the finite differences by derivatives. The 
result will read 

{50(A), 50(m)} = 0, 
{5+(A),5+(m)} + ^+(A)5+(/x)-5+(/z)5+(A) = 0, 
{S^iX),S^if,)}-Si{X)Srif^)-Srif^)Si{X) = 0, 

{s^{x),snf^)} + s?Wst{f^) + ^^^^^ = 0, 

{S?(A),5r(/i)} + 5?(A)5r(M)-^^^4^^ = 0, 

A — jJi 

{5r(A),5+(/x)} + 5r(A)S+(/i)-5+(/.)5r(A) + 2^^^^4^-^ = 0, (22.6) 

where the dot denotes the derivative with respect to t. We call an algebra defined by these 
relations the classical quasi-Gaudin algebra. The existence of this algebra can easily be proved 
if one replaces n by nTi in formulas ( |5.2| ), takes the limits 7i — > and n — > oo, and checks that 
the resulting generators (having the same form as in ( ^.2| ) but with n replaced by t) obey the 
relations (22.6). Taking the same limits in formula ( |22.4| ), we obtain 

H,iX) = 5°(A)50(A) - 5+(A)Sr(A). (22.7) 
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It is not difficult to see that functions ( 22. 7| ) are in involution for any given t, 



{Ht{X),Htif,)} = 0, (22.8) 

but this is generally not true if the values of t differ from each other. This means that, exactly 
as in the quantum case, the classical quasi-Gaudin algebra generates an infinite number (or, 
more exactly, a continuous family) of integrable classical systems. It is natural to call them 
the classical quasi-Gaudin models. 

According to the famous Liouville theorem, all these models are integrable in quadratures 
and, in this sense, can be qualified as exactly solvable. Remember, however, that the practical 
usefulness of the Liouville theorem is very low if one deals with classical systems with many 
degrees of freedom. In general, we simply cannot present their explicit solutions!|^ At present 
time there is only one constructive way for solving classically integrable systems. It is given 
by the famous classical inverse scattering method (CISM) (see e.g. [p]]) which presents more 
or less explicit algorithm of constructing action-angle variables for models associated with 
classical Yang-Baxter algebras]^. We expect that the classical quasi-Gaudin models which 
we obtained are in fact quasi-exactly solvable in the sense of CISM. This means that using 
the standard prescriptions of this method and applying them to models ( 22.7] ) we probably 



will obtain only some part of action-angle variables, but not a complete set of them. It 
would be very desirable to try to verify this conjecture because if it is true then one obtains 
obtain a remarkable possibility of introducing a concept of quasi-exact solvability in classical 
mechanics. 

Note that the relations (|22.6| ) can be rewritten in more standard form if one introduces a 
new bracket: {{ , }}. Defining it as 

{{at, bt}} = {at, bt} + g{at)htht - g{bt)btat, (22.9) 

where g{a), as usually, denotes a grading, we obtain for ( p2.6|) : 

{{5°(A), 5° (/.)}} = 0, 

{{S+iX),S+ip)}} = 0, 

{{S-iX),S-{f,)}} = 0, 

{{50(A), 5+ (/x)}} = -(A-/x)-H5+(A)-5+(^)}, 

{{5°(A), 5- (/.)}} = +(A-/.)-H5-(A)- 5- (/.)}, 

{{5- (A), 5+ (/.)}} = -2(A-m)-HS°(A)-50(^)}. (22.10) 

We see that the classical quasi-Gaudin algebra can be considered as an ordinary classical 



Gaudin algebra realized with respect to the non-standard bracket ( 22.9 ). We arrived at the 



same situation which we had in section 15 for quasi-commutators. This similarity is not 



^*This situation is exactly the same as in quantum mechanics: we become happy when a certain quantum 
problem is reduced to a purely algebraic one and call such a problem exactly solvable (or quasi-exactly 
solvable). But it does not mean that we are able to solve any algebraic problem explicitly, especially if the 
number of unknowns is large. 

^^The fact that CISM works for integrable models associated with classical Yang-Baxter algebras does not 
mean that it works for any integrable model. 
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accidental because the quasi- Poisson bracket ( 22. 9| ) is nothing else than the classical limit of 
a quasi-commutator which can be defined as 



[[an-, K]] = an+g{h)K - K+g{a)0'n- (22.11) 

Indeed, following the standard prescriptions of dequantization procedure, we should rescale 



the relation (22.11) as 

[[anh, bnh]] = 0-nh+g{h)hbnh — bnn+g{a)hO'nh (22.12) 
and then rewrite it in the form 

[[ann, bnn]] = Ka, b^n] + h^^^^±^^f^^b^n - n^-^^^±^^f^a^n (22.13) 

Taking now the limits 7i — > and n ^ oo, introducing a new finite and continuous variable 
t = nh and using definition ( 22. 9[) , we obtain 



[[ann,bnn]]^H{at,bt}}. (22.14) 

From this analysis it follows that the definition ( ^2.91) works not only for the objects with 
gradings and ±1 but for arbitrarily graded ones. One can check that the quasi- Poisson 
bracket is anti-symmetric 

{{at,bt}} + {{bt,at}} = (22.15) 

and obeys the standard Jacoby identity: 

{{ {{at, h}}, ct}} + {{ {{bt, q}}, at}} + {{ {{ct,at}}, h}} = 0. (22.16) 

For zero- graded objects it coincides with an ordinary Poisson bracket 

{{at,bt}} = {bt,at}, if 5(a) = 5(6) = (22.17) 

and therefore the commutativity of functions Ht{X) can be considered as an elementary 
consequene of the commutativity of their Gaudin analogues H(X). 

23 The perspectives 

We completed the exposition of the proposed i?-matrix approach to the problem of quasi- 
exact solvability. Summarizing, we can say that we found a simple modification of the 
celebrated quantum inverse scattering method which leads to new wide classes of completely 
integrable models. The latter are distinguished by the fact they admit an algebraic Bethe 
ansatz solution only for some limited parts of the spectrum, but not for the whole spectrum, 
as in the standard version of QISM. Therefore they are typical quasi-exactly solvable models 



discussed in refs. pi| , 45]. An underlying algebra responsible for both the phenomena of 
complete integrability and quasi-exact solvability was constructed. We called it the "quasi- 
Tr algebra" and demonstrated that it can be considered as a very simple (but non-trivial) 
deformation of the ordinary Tr algebra of monodromy matrices generated by various solutions 
of Yang-Baxter equation. This enabled one to claim that we actually have found a new class 
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of fundamental symmetries of integrable quantum systems. The main attention in the paper 
was devoted to construction of such deformations for simplest Tr algebras, i.e. for Gaudin 
algebra ^[s/(2)] and its quantized version — the Yangian 3^[s/(2)]. We also constructed and 
solved the quasi-exactly solvable models of magnetic chains generated by these algebras in 
the framework of QISM. The classical vesions of quasi- Tr algebras were discussed briefly and 
we intend to return to this interesting question in the nearest future. 

The main result of this paper lies however in the proof of the fact that any graded asso- 
ciative algebra with unit admits such "quasi" deformation. We showed that this deformation 
modifies the structure of the elements of an algebra and changes the definition of their asso- 
ciative product. At the same time, despite the numerous attempts we did not manage to find 
any appropriate deformation of a co-product which would preserve a possible Hopf algebraic 
structure of the initial algebra. This enables one to make a conjecture that quasi- Tr algebras 
which we constructed in this paper are not the Hopf algebras^ 

The method we proposed in this paper is interesting primarily because it opens a practical 
way of constructing quasi-exactly solvable deformations of all known exactly solvable and 
completely integrable quantum models obtainable in the framework of standard version of 
QISM. The list of such models includes: 

1) Exactly solvable models of ordinary quantum mechanics. In one-dimensional case 
these are: the simple harmonic oscillator, the Kratzer and Morse potentials, Poschel-Teller 
potential wells, etc. In multi-dimensional case these are: the multi-dimensional harmonic 
oscillator and also various separable models defined on non-trivial curved manifolds and 



constructed in ref. [^]. As was demonstrated in [46|, all these models can be considered as 
various degenerate cases of Gaudin models and hence can be obtained in the framework of 
classical i?-matrix approach. This, in turn, means that we can apply to these models our 
deformation procedure and obtain their quasi-exactly solvable conterparts. This would be 
a perfect test for the proposed method because the explicit form of all these quasi-exactly 
solvable models is already known (see e.g. refs. ||4^, p6t|). 

2) Exactly solvable models of multi-particle quantum mechanics. These are: the Calogero 



model, the Calogero-Moser model and other similar models discussed for example in ref. |24]. 
Recently it was demonstrated that these models can be obtained in the framework of the so- 
called method of dynamical i?-matrices. A main distinguished feature of this method lies in 
the fact that the i?-matrices used in it depend themselves on dynamical variables, i.e. on 
the components of a monodromy matrix. For us however it does not matter what kind of 
i?-matrices is used in commutation relations of algebra Tr because of the universality of the 
proposed deformation method. This means that we can use this method for constructing 
quasi-exactly solvable versions of multi-particle Calogero and Calogero-Moser models. Note 



that particular classes of such models have already been constructed in ref. |44] in different 
way (see also the book [^]). 

3) Exactly solvable models of spin chains. These include, first of all, the local chains: 
the simples 1/2-spin chains with XXX, XXY and XYZ symmetries as well as their higher 



spin generalizations constructed in refs. |35, 34 1, and also the various inhomogeneous spin 



chains with the same symmetries. These models are associated respectively with rational. 



^"it would be nice to call such algebras "quasi-Hopf algebras", but, unfortunately, this term is already used 
in a different sense 
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trigonometric and elliptic solutions of Yang-Baxter equation. In this paper we discussed only 
a very particular case of XXX models. Moreover, for these models we presented only the 
sets of commuting integrals of motion but did not discuss the methods of extracting from 
these sets the hamiltonians of physically interesting systems. 

4) Exactly solvable models of field theory. We know a lot of models of such a sort 
including non-relativistic models like non-linear Schrodinger equation, sine-Gordon model, 
various lattice models, and also the relativistic models describing various forms of four-fermion 
interaction. Using our deformation technique, one can try to construct their quasi-exactly 
solvable conterparts. This, may be, requires some comments. The quasi-exactly solvable 
problems which were considered up to now were defined as problems having a certain finite 
parts of the spectrum constructable in an algebraic way. If the number of such algebraically 
constructable states is A^, then, in order to find them, one should solve a certain A^-th order 
algebraic equation. A fruitfulness of such a concept in quantum mechanics is very well known. 
However, if we want to have something like field-theoretical quasi-exactly solvable models (or, 
simply, quasi-exactly solvable models with many degrees of freedom) one should allow A^ to 
tend to infinity, because the only physically interesting situations in such models may appear 
when the number of excitations is infinitely (or, simply, very) large. But for this we should be 
able to solve the corresponding "algebraic" equations! If we speak only of the algebraizability 
itself, as, for example, in the partial algebraization method, then we are unable to solve this 
problem, since the complexity of corresponding algebraic problems catastrophically increases 
with the increase of A^. However, if we have some additional structure in our systems, which 
is, say, the integrability property, then, as we have seen above, the equations for solving our 
algebraic problem take the form of Bethe ansatz equations which in the large A^ limit are 
well defined and have been studied in detail in many occasions. 

5) Exactly solvable models of statistical mechanics. We know that the mathematical 
structure lying behind such models is again the Tr algebra, which however has now the 
meaning of an algebra of transfer matrices. It would be very interesting to try to understand 
what kind of statistical systems correspond to quasi- algebras, and in which sense they are 
quasi-exactly solvable. 

Of course, the construction of quasi-exactly solvable models with an infinite number of 
degrees of freedom in the framevork of quasi-QISM approach is a long-term program which 
hardly can be solved in the nearest future. Now we can only express an optimistic hope 
that its realization is possible. The reason for our optimism is based on two facts: a) that 
the standard QISM formalism is well adopted to the study of exactly solvable models with 
many degrees of freedom, and b) that there is no essential difference between QISM and 
quasi-QISM, as it follows from the results of our paper. Of course, the quasi-QISM does not 
exactly coincide with the standard QISM and has many specific features which need a very 
careful study and understanding. In this connection, I would like to mention two important 
problems which immediately arise if one wants to give a physical interpretation of the results 
of the present paper. 

The hermiticity problem. Any physically realistic model should be described by her- 



mitian hamiltonians. If one looks at quasi-Gaudin and quasi-XATX models (6.1) and (12.1) 
which we constructed in our paper, then it becomes obvious that their hermiticity properties 
strongly depend on the choice of a representation in which the corresponding ^[si(2)] and 
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3'[s/(2)] algebras act. One of trivial possibilities of making these hamiltonians hermitian 
and the corresponding wavefunctions normalizable is to take the unitary representations of 
these algebras which, due to their incompactness, are all infinite-dimensional. Are there any 
other ways? The experience accumulated in the study of integrable and quasi-exactly solv- 
able models enables one to answer this question positively. Indeed, there is a constructive 
method for building such models (and this is just the inverse method of separation of vari- 
ables discussed in section ^) which allows one to guarantee the hermiticity of the resulting 
hamiltonians from very beginning, i.e. by construction (for more detail see e.g. ref. p5| , |4^). 



Unfortunately, in order to apply this method to models ( |6.1| ) and (12.1) as well as to other 
models obtainable in the framework of quasi-QISM, one should first separate the variables in 
them (we mean here the separation of variables in the generalized Sklyanin sense ||3^ ) . This 
immediately creates two important questions: a) Is there any regular way for performing 
Sklyanin's procedure of separation of variables for models obtainable in the framework of 
quasi-QISM? b) Is there any direct way of formulating a criterion of hermiticity immediately 
in terms of quasi- Tr algebras and their representations? The second question creates, in turn, 
a problem of developing some kind of representation theory for quasi- Tr algebras, etc. 

The problem of locality. This is no less important problem on whose positive solution 
will finally depend a progress in the realization of our program. We know that most of 
physically interesting models of field theory are local. At the same time the models (|6.l| ) and 



(12.1) look as typically non-local ones. Then the following question immediately arises: Is it 
possible to extract from the set of integrals of motion of these models some local hamiltonians? 
And, if not, then may be there exist other realizations of quasi-^[sZ(2)] and quasi-3^[sZ(2)] 
algebras for which it would be possible? 

I hope that the hermiticity and locality problems will finally find their solutions in the 
framework of quasi-QISM approach. The more that the construction of hermitian and local 
quasi-exactly solvable models with many degrees of freedom is possible even in the standard 
version of QISM. I am addressing the interested reader to ref. in which this question is 
discussed in detail. 
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